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ABSTRACT 

The use of hydraulic fracturing in the field of energy 
resources development has generated a aaed for effective 
methods to determine the size and orientation of induced 
fractures. Surface deformation monitoring has been 
considered as one of the feasible geophysical methods. In 
order to understand its relevance, we have to understand how 
well the surface deformation delineates changes of the 
source parameters of fracture. 

I review past studies of hydraulic fracturing; field 
experiments of deformation monitoring and basic problems in 
the behavior of hydraulic fractures. Mainly using an 
analytical solution for the surface displacement due to a 
horizontal penny-shaped fracture, I discuss a method to 
determine the source parameters from surface deformation 
measurements, and present a measure to evaluate the 
applicability of Rootes renon it oriia for investigation of 


hydraulic fractures. 
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1. Introduction 

Since its introduction to the petroleum industry (Clark, 
1949), the technique of hydraulic fracturing has _ been 
applied as a tool for underground waste disposal, geothermal 
energy recovery, Solution mining and in-situ bitumen 
recovery from oil sands and oil shale deposits as well as 
for conventional oil/gas reservoir Stimulation. The 
technique has received a considerable amount of theoretical 
and experimental attention. Unfortunately, the analysis of 
the behavior of induced fractures is still inadequate. The 
difficulty is that the form of the fracture depends on 
several variables such as the interaction between the 
properties of injection fluid and the properties of the 
formation rock, the strength characteristics of the 
reservoir materials, the in-situ stress state, the initial 
pore pressure and the fluid penetration, and the 
pre-existing fractures and zones of weakness. 

Hydraulic fracturing in a project should be designed to 
create fractures of a desired form for the particular 
purpose. The form of the Fracture, especially its 
Orientation and size, is crucial. Nevertheless there have 
been few cases in which the form of induced fractures could 
be well determined. One of the biggest problems in the 
application of hydraulic fracturing relates to the methods 
for the determination of the size and orientation of 
fractures. Geophysical techniques involving surface 


measurement of physical quantities and their changes 





noispubewsynl .t ; : 


iefD}. etsevinl.  muefentges. sha ov neresubougse est eonte 

















host 1uszs3? siluszbyd «Fo -supindses, edz. (TERT 


set Cesare th pat STEN YE 4 isc3 Bs > 
so roe tbh -etasv. bAdarzeyghnc 16h ioc? & 28 bs i iggs 7 
upgie-nt ~ dhs @ SAinin Mocs .ee ,vieveos1 — yp Tene 


j ra -* 
az >4i:coueb.elsaz! Liotbngs shasz .f¢ mors y7e 


LSciseues “iiewé sidset#iteres « bavieos: zed supindses 
0 
ey lions sa plesehuicsIein,nowtnesis Jaetnemizeqxe @ 
Vs «<a 

fe @! eSimwoessl Besuber to tolvaded ae 


7 


12, Sav sae viet eae of jed3 ei yiluoliaie 


ijnasvsetnt Sia 2a foue asidaiisvy isis 


2 2aI3e US271 2.4/135768 3i0TI8 
$185 aed eh, an? ons S IVE039 


2@enisew 3© eens: bone eaqweset? paisetxec 


29f ot Sluene toefeda & at enlavsogs2 ai ivszvavis 


ry 
w 
~~ 
W 
F 
t 
ai] 
<4 
iB 
‘i 
~ 
Cc 
‘ys 
ua 
oD 
of 
=) 


} 
ge) ica ha et @ (aie oe on Al its : ~ 
edi vilseioeges Siudaaoa sis 20.7023. e017, ,sa0qu7 


és 
. 
rt) 
i? 
tg 
rat 
Lee 


hs “ PS e ’ iy 5 ‘\ = = 
a@eeisthsisvet .isicuts 2 ,ente Sas noitstnel® 





: 
- a 
- 


vine re ee 


associated with hydraulic fracturing are required. Any 
technique sensitive to the presence of fractures should be 
Sea acedh The feasibility of seismic and acoustic 
measurements, ground deformation monitoring and electric 
potential measurements have _ been evaluated by field 
experiments, but the number of the experiments with 
published results is quite small. Here I discuss the use of 


ground deformation measurements. 


1.1 Measurement and Analysis of Ground Deformation 

The value of ground surface deformation measurements 
for the determination of quantitative details of fracture 
implies the detectability of small magnitude surface 
deformation and methods for the analysis of the resulting 
data. Accurate measurements with high resolution instruments 
and construction of. a tractable, realistic fracture model 
are the minimum requirements. 

Even if surface deformations are measured accurately by 
high resolution instruments, this does not imply the 
detection of the surface deformations due to the fluid 
injection into the subsurface. The earth is not a stationary 
but a deformable body which changes in internal structure 
and shape by natural and induced forces. The types of 
surface deformation are multifarious in terms of time scale, 
Spatial extent and magnitude. The extraction of the signals 


from the fluid injection operation can not be attained 
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without suppression and estimation of the background noise 
and proper data processing for the signal enhancement. 

The major noise which should be taken into 
consideration for the extraction of the signals from the 
fluid injection is caused by tides and meteorological 
effects such as diurnal thermal strain, wind, barometric 
fluctuation and moisture change, and the near-surface 
enviroments such as the machinery for injection operation 
and fluctuating aquifers. Only some of them are periodic and 
predictable. (Wyatt and Berger, 1980; Savage and Prescott, 
1973). 

Hydraulic fracturing is not always accompanied by 
permanent deformation of the ground surface. The period of 
the fluid injection operations usually ranges from _ several 
tens a minutes to a few hours, or a few days. The 
measurements are required to follow relatively short time 
phenomena. 

In order to analyze the observed surface deformation, 
we assume a mathematical fracture model, and we need to 
understand the theoretical ground deformation due to the 
model fracture. There are two ways to characterize the 
subsurface deformation. In one of them the subsurface 
deformation is characterized by. the displacement 
dislocation, and the other assumes the forces acting in the 
region of subsurface deformation. The complexity of real 
subsurface deformations makes it difficult to suggest how 


the distribution of the displacement dislocation or the 
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forces should be modelled. Surface deformation associated 
with fluid injection has not yet been studied sufficiently 
that a strict fracture model can be constructed. 

We need a tractable and realistic crack model in order 
to start the solution of this problem. Even if the crack is 
assumed to have simple shape characterized by uniform forces 
Or constant displacement dislocation in an idealized elastic 
earth model, analytical solutions exist only for a limited 
number of symmetrical geometries ( e.g., Sih,1973; Kassir 
and Sih,1975 ), and even those are complicated. These 
solutions all have a singularity, a stress concentration of 
Paeanity, ac the crack tip, which probably plays an 
important role in the propagation of the fracture but is 
physically impossible. Therefore the assumption of the 
displacement or stress conditions should be reconsidered to 
make the crack model more reasonable. For more complex crack 
models, numerical approaches are naturally required. It is 
also important to understand how much the _ theoretical 
surface deformation fields from the simple crack model 
differ from the alternative model which takes into 
consideration ‘such a physically impossible situation, and 
how well the surface deformation delineates the changes of 
the source parameters of fracture. 

There have been few hydraulic fracturing experiments 
where ground deformation was measured and the results 
published. The reason for this seems to be the difficulty of 


constructing a tractable, realistic crack model and of 
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developing the techniques for field measurements. 


1.2 Source Models in Volcanism and Seismology 

There have been several source models applied to the 
analysis of surface deformation associated with volcanism. 
Yamakawa (1955) and Fukuda (1944) assumed spherical models. 
Walsh et al. (1971) studied a semi-infinite line source 
model. Davis et al. (1974) used an ellipsoidal, anelastic 
source model for the interpretation of the Riiaces (s 
observation. Dietrich et al. (1975) calculated the surface 
displacements due to a variety of volcanic source shapes by 
finite elements analysis. Swanson et al. (1976) applied the 
above five different models to the observations around 
Kilauea's rift and compared the results. Koide et al. (1975) 
discussed formation of fractures around magmatic intrusions. 
Lachenbruch (1961) analyzed natural extension fractures 
which initiate at the surface and propagate downward. 

‘In seismology, elastic theory of dislocation 
demonstrated the equivalence of small shear dislocation and 
the double-couple force, the SO called body force 
equivalence (Steketee,1958; Maruyama,1963; Burridge and 
Knopoff, 1964), and then brought a significant advance with 
regard to the elucidation of seismic source mechanism. The 
relation of seismic source to shear fault movement has been 
substantiated by many observations. Many earthquakes have 


been described in terms of the time and spatial parameters 
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of such seismic fault models, (cf., Yonekura, 1978). It is 
furthermore noted that surface deformations and seismic 
radiation have been explained quite reasonably by this 
rather simple model. 

Basic mathematical treatments of elastic theory of 
dislocation for the calculation of surface deformation were 
presented by Steketee (1958) and Maruyama (1964). Analytical 
solutions for simple shear fault models with uniform 
displacement dislocation on a plane were derived by Chinnery 
(1960), who assumed a two-dimensional, vertical, strike-slip 
fault, Press (1965), who analyzed a vertical, rectangular 
fault model with dip and strike-slips (and published 
incorrect results), and Manshinha and Smylie (1971), who 
took account of arbitrary dip-angle of a rectangular fault 
plane. The static deformation field due to dislocation in a 
multi-layered medium has been studied by Sato (1971) and 
Serommeande Matsuurs (1973.1..1974,., 1975). Nyland (1971) 
discussed body force equivalents as sources of anelastic 
processes. Other references are found in articles of 
Ben-Menahem and Singh (1969, 1970). 

Several attempts have been made to interpret static 
displacement fields due to the occurrence of an earthquake 
in terms of a dislocation function varying along a complex 
fault surface. In those analyses, the dislocation function 
has been determined by using the method of linear inversion, 
the least-squares procedure or the finite elements 


modelling, under an assumption that the other presupposed 
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fault parameters are adequate. The reliability of the 
estimated results depends strongly onspthat,+,of:.the 
presupposed fault parameters. (Savage and Hastie, 1966; 
Canitez and Toksoz, 1972; Mikumo, 1973; Jungles and 
Franzier,1973; Alewine and Jordan,1973; Nyland,1973, 1976; 
McCowan et al.,1977; Miyashita and Matsuura,1978). 

Earthquakes represent failure of geological materials 
to resist tectonic stress. The stress is accumulated over a 
period of time, and must be accompanied by deformation and 
Straining of rock. Laboratory experiments and theoretical 
considerations predict that a change in the strain rate 
should precede the failure. Therefore, it is a crucial 
subject for earthquake prediction to monitor the temporal 
and spatial extent as well as the magnitude of the changes 
in the strain rates preceding earthquakes. 

The value of the deformation measurement is also found 
in the fact that a substantial part of the faulting 
associated with earthquakes takes place too slowly to be 
detected by most seismometers. There is | the further 
possibility that slow and 'silent' earthquakes may result in 
Significant displacements, and anelastic processes following 
earthquakes may strongly affect the stress enviroment. These 
phenomena would be accompanied by episodic ground 
displacements that might not be much smaller than those 
associated with seismic events. Geodetic monitoring of 
tectonic deformation has been conducted in various forms, 


for example, the monitoring of long-term strain accumulation 
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patterns and rates, medium-term episodic strain changes as 
the preparatory phase of an earthquake, the detection of 
short-term strain changes just prior to its fault pease rial, 
and the observation of postseismic deformation as _ the 
relaxation process of regional and local stress patterns. 
Not only the -measurement of local strain accumulation and 
release in and around seismic active zones but also global 
tectonic deformation monitoring play a great role in the 
study of the dynamics of the earth interior. ( NRC, 1981; 
metanc, tatvey Thatcher, Th9795 ) Savage sceteral #71981; “J.G.R 
Vol.85,No.B4,1980). 

The interaction of a fracture and a free surface has 
been studied also in the field of fracture mechanics; e.g., 
Sih (1973), Erdorgan et al.(1974). Stress intensity factors 
on the fracture periphery are of much interest in fracture 
mechanics, and stresses and displacements on the free 


Surface are not given generally. 


1.3 Previous Field Work 
I refer to four field experiments; (1) Ren Jen Sun, 
1969, (2) Holzhausen et al.,1980, (3) Evans et al.,1980 and 
(4) M.D.Wood. Inc.,1979, in which the behavior of hydraulic 
fractures was analyzed by using the surface deformations 
measured by either levelling or borehole tiltmeters. 
Previous field experiments suggest that fluid injection 


into the subsurface certainly causes surface deformations, 
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and furthermore demonstrate that the measurement of surface 
deformations can provide an indication of the overall 
geometry of hydraulic fracture and reflect major changes in 
the facture behavior such as a change in propagation 
direction (e.g., from vertical to horizontal) or cessation 
of growth resulting from intersection of pre-existing 
natural joint or fissure. 

The magnitude of the measured surface deformations is 
large enough compared with the resolution of the instruments 
and the survey accuracy. The tiltmeter surveys reported that 
tilt changes of ground surface are on the order of 0.1-2.5 
microradians, (cf., (2), (3) and (4) ), whereas the 
resolution of the tiltmeters used for the field experiments 
is order of 10°? - 10°? microradians. The uplift survey in 
one of the field experiments reported maximum surface 
wees OL) 4-12 mm. 4s cEe, (Ne JRO wUnfortunately, the 
specification of the levelling surveys is not clear. 
Probably the best attainable measurement accuracy OL - 
levelling survey is estimated by 4mm x (line length in km)°°, 

The period of the fluid injection operations ranged 
from several tens of minutes to several hours, and for one 
case, (2), it was as long as about 40 days. Even if a 
fracturing operation was conducted over a long term, ground 
responses were reported to be characterized by several 
short-term events which suggest the repetition of growth and 
collapse of peertec tires: ines Grea’ size of surface 


deformation around the well is estimated to be roughly 2-3 
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times of the injection depth (180-350 m). As a reference, 
the calculated size of the fractures ranges from one-fourth 
to half of the injection depth in most cases. Although 
results have not published, the tiltmeter surveys are 
reported to have been conducted for hydraulic fracturing as 
deep as 3000m. (M.D.Wood.Inc., personal communication). 

Surface deformation in the previous experiments was 
measured by levelling or by tiltmeters. The previous tilt 
Surveys consisted of at most nine borehole tiltmeters, eight 
of which were installed in a circular array from the well at 
a radius of half the injection depth. The number of the 
tiltmeters used does not seem to be sufficient. 

In order to measure complex surface deformation more 
exactly and delineate the fracture behavior more clearly, we 
need to incorporate more observations of the measurable 
quantities. On the other hand, we have to examine the 
techniques of the measurement from the aspect of the 
relation between the the accuracy of each survey and the 
predicted magnitude of surface deformation, the measuring 
efficiency and the sensitivity of each measurable quantity 
to the form of the fractures. We could use optical or laser 
levelling, Surface laser ranging and spaceborne laser 
ranging for the displacement measurements. The disadvantage 
of these surveys is the difficulty of continuous operation. 
For strain measurement, several strain meters are available, 
such as long baseline Invar wire or fused quartz rod strain 


meter. For tilt measurement, several types of tiltmeters are 
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available. In addition to tilts and uplifts, measurements of 
horizontal displacements and strains might be useful. and 
should be considered. 

The mathematical fracture model so far used for 
analysing the observed surface deformations is either a 
-horizontal penny-shaped model, (Sun, 1969), Or a 
two-dimensional inclined model, (Pollard et al., 1979). Both 
models assume a plane fracture subjected to a uniform 
internal pressure in an idealized earth such as a 
semi-infinite, homogeneous and isotropic, linear elastic 
medium. In spite of the simplification, analyses of the 
observed surface deformation led to the plausible 
conclusions about the fracture behavior. The two-dimensional 
inclined fracture model was applied to one fracture among 
Sixteen fractures whose corresponding surface deformations 
were analysed quantitatively in the previous field 
experiments. 

In one of the field experiments, (Sun, 1969), the 
radius of the fractures was estimated from the volume of the 
injection fluid. For such an estimate of fracture radius, 
the values of the elastic constants and the internal 
pressure should be known, and an assumption about "che 
impermeability of surrounding rock should hold true. Such 
analysis is not always applicable to all cases. The 
comparison between the observed and calculated uplifts 
Supports the validity of the mathematical fracture model. 


The calculation of the fracture size from the observed 
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uplift data, so called "inversion problem", was not 
discussed at all in this project. 

In the other three field experiments with tiltmeter 
measurements, the size of the fractures was estimated from 
the observed tilt fields by two methods. The theoretical 
surface deformation field based on the fracture models is 
expressed by the product of the quantity given by elastic 
constants and internal pressure in the fracture and the 
quantity given by functional See between the 
measurement point and the source geometry. The former 
Quantity will be called here "the amplitude factor” of the 
surface deformation. If the geometric relation between the 
measurement points and the source allows, we can estimate 
both of the amplitude factor and the geometric parameters of 
the fracture from the observed surface responses. 

When the induced fracture is regarded as horizontal and 
ea honeds eight tiltmeters installed in a circular array 
around the injection well, ( such site configuration of the 
instruments has been used in most previous experiments), 
provide substantially only the tilt responses at one 
meaSurement point. In such case, the amplitude factor should 
be known a priori to permit estimation of source parameters. 
Nevertheless, the physical quantities which contribute to 
the amplitude factor are very uncertain. An estimate of the 
internal pressure in the fracture (driving pressure) is very 
difficult. The uncertainty of the driving pressure has been 


overcome by using the concept of stress intensity factor. 
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Linear elastic fracture mechanics states that a fracture 
propagates when the stress intensity factor at the fracture 
tip reaches a critical value (called fracture toughness) of 
the material. Then we have a relationship between the 
driving pressure, and radius of the fracture and the 
fracture toughness of the surrounding rock. It is also a 
fact that estimation of the driving pressure from the 
relationship requires a certain value of the fracture 
toughness and the radius, which are usually unknown and 
Should be calculated from the observations. However, 
assumptions about the order of the magnitude of fracture 
Size and the fracture toughness gives a rough estimate of 
the driving pressure. Once the amplitude factor is 
estimated, it is easy to determine the fracture radius. 
Except for one example in which the observed wellhead 
pressures were used for the analysis, there is no 
application of the wellhead pressure to estimate the driving 
pressure. 

The other method (Evans et al.,1980) for determining 
the geometric parameters of a fracture was based on the 
two-dimensional inclined fracture model. The eight 
tiltmeters installed in a renee array provided 
information to determine both of the amplitude factor and 
the geometric parameters such as the strike, dip-angle and 
length of the fracture. If the observed tilt field is 
normalized by the magnitude of a reference measurement 


point, the uncertainty of the amplitude factor is estimated 
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from the absolute magnitude of the measured surface 
deformations. Once the geometric parameters of the fracture 
are determined from the relative magnitude of observed 
tilts, they lead to the amplitude factor. The driving 
pressure estimated this way is valuable for understanding 
the pressure distribution in as fracture and the in-situ 
fracture toughness. This analysis is essential when compared 
with the previous one which depends heavily on uncertain 
assumptions. | 

In the tiltmeter surveys, the surface tilt responses 
have been monitored relatively densely in time. If the 
surface deformations had been analyzed in terms of the time 
history of the radius and driving pressure of the fractures, 
Significant results which delineate the process of fracture 
growth might have been obtained. As described in Section 
mete tneorles Of tractUre propagation predict ~ that a 
fracture is not necessarily a plane fracture centered at the 
injection depth. It seems that, in the previous experiments, 
the data acquired under the limitation of the small number 
of instruments and the site configuration did not stimulate 
the investigators to verify the possibility of those 
theoretical predictions. 

Except ror the use of observation wells, most 
diagnostic tecnologies for the inspection of induced 
fractures, such as wellbore impression packers, injection of 
radioactive sand into the fracture for subsequent detection 


with a well logging tool, lineament analysis, seisviewer 
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imagery, borehole television and post-fracture well 
temperature logs, suffer the disadvantage that they give, at 
best, Practere dimensions and orientation at the wellbore 
only. The fact that past field experiments of surface 
deformation measurements, somehow or other, delineated the 


overall size of the fractures is valuable. 


1.4 Plan of the Thesis 

In this thesis, I first review the four field 
experiments cited above, and some _ basic problems Ce 
hydraulic fracturing; the initiation and orientation changes 
of fracture. Next, I discuss theoretical size of surface 
deformations, using a horizontal penny-shaped fracture 
model, inclined tensile and shear fracture models. Next, an 
inversion scheme relating the surface deformations due to a 
horizontal penny-shaped fracture is introduced. The 
resolution inherent of this inversion scheme is evaluated 
against an eStimated level of measurement errors. The 
optimum survey configuration and a possible combination of 
different observable data are discussed from the point of 
view of the attainable highest resolution. Assuming that the 
surface deformations during an injection operation are given 
by a time function, an effective procedure for the 
determination of rates of fracture growth i$ shown. 
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2. The Behavior of Hydraulic Fracture 
First, I refer to four previous field experiments; Ren Jen 
Sun (1969), Holzhausen et al. (1980), Evans et al. (1980) 
and M.D.Wood. Inc. (1979), in which surface deformation 
measurements were used to investigate hydraulic fractures. 
The major content of these field experiments is tabulated in 
Table 2.1.1 ( at the end of Section 2.1 ). In the next two 
sections, the basic problems of hydraulic fracturing; the 
initiation and orientation changes of a fracture, are 


reviewed. 


2.1 Previous Field Experiments Utilizing Deformation 
Measurements 

(1) Ren Jen Sun (1969) 

Hydraulic fracturing experiments were conducted at ORNL 
(Oak Ridge National Laboratory, Tennessee) from 1960 to 1965 
to evaluate the feasibility of disposal of radioactive waste 
by injecting grout mixtures into a shale formation. Through 
the experiments, the theoretical size of hydraulically 
serendé fractures was compared with field measurements of 
the surface uplifts. ORNL is located in the Appalachian belt 
of faulted and folded Paleozoic rocks. All the injections 
were made into the lowermost unit of the Conasauga shale 
formation, which is about 300 m thick. 

The experiments consisted of seven injection operations 


and five uplift surveys. The behavior of one of the induced 
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fractures is as follows. The injection was made through a 
Slot in the casing at a depth of 285 m. The estimated 
overburden pressure at the injection depth is 7.57 MPa. 
Fracturing was initiated with water. The fracturing started 
at 10.39 MPa of pressure at the surface, then the pressure 
fell quickly to 9.31 MPa at a pumping rate of 0.57 m?’/min. 
and increased to 11.27 MPa as the rate was increased 
slightly. Six minutes after the fracture was started, by 
which time 3.2-3.4 m* of water had been pumped, the pressure 
in a nearby observation well, 9 m west of the injection 
werl,) suddenly rose to 5.49 MPa, indicating that the 
fracture had extended into it. After the pump was. stopped, 
the water was bled back out of the injection well, and the 
pressure in the observation weli fell to 2.45 MPa in about 
100 minutes. The main injection of water-cement-bentonite 
mixture was then started. The mixture consisted of 346 m°* of 
grout containing 201 tons of cement and 5700 kg of bentonite 
clay. The maximum pressure for the slurry injection was 
15.88 MPa. About 40 minutes after the injection started, the 
pressure dropped to between 11.76 and 12.05 MPa. 

The history of pressure during the injection is shown 
in Fig.2.1.1. The fact that the injection pressure was, in 
all cases, greater than the formation overburden pressure, 
indicates the occurence of horizontal fracture. Study of 
cores made after grout injections also showed that the 
fractures were concordant with the nearly horizontal 


bedding. Fig.2.1.2 shows the calculated and observed surface 
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PUMPING RATE. IN CUBIC METERS PER MINUTE 
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history during a _ hydraulic fracturing 
operation (from Ren Jen Sun, 1969). 





: EXPLANATION j 
@ Surveyed uplift along line between benchmarks 4A ond 100 
© Surveyed uplift along line between benchmarks 3C and 3E 
4 x Surveyed uplift along line between benchmarks 68 and OF 
x 
0 : 
700 400 300 200 100 100 200 300 400 500 600 700 








RADIAL DISTANCE BETWEEN SURVEY BENCHMARKS AND THE INJECTION WELL, IN METERS 


Fig. 7. Calculated and surveyed uplift produced by injection 1. second experiment. Septem- 
ber 3, 1960. 


EXPLANATION 
& Surveyed uplift along line between 
benchmarks 4A and 100 
© Surveyed uplift along line between 
benchmarks 3C and 3€ 


* Surveyed uplift along line between 


400 300 200 100 100 200 300 400 500 600 700 
RADIAL DISTANCE BETWEEN SURVEY BENCHMARKS AND THE INJECTION WELL, IN METERS 


Fig. 8. Calculated and surveyed uplift produced by injection 2, second experiment, Septein- 
ber 10, 1960. 


Piged. Ices ate samp les Ofoica lomlated and 
surveyed upli;is (from Ren Jen Sun, 1969). 
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uplifts. Fig.2.1.3 shows the observation geometry, the 
measured surface uplift and the extent and thickness of the 
injected grout sheet obtained from core-hole data for two 
cases of the experiment. 

The calculation of the theoretical uplift is based on a 
mathematical model of a horizontal penny-shaped fracture 
internally pressurized by fluid injection. For the model 
fracture in an infinite isotropic and homogeneous, 
impermeable, linear elastic medium, the fracture radius, a, 
can be calculated from the internal pressure P, the total 


volume of injection fluid Q and elastic constants, say 





Young's modulus E and Poisson's ratio P , by the relation 
SEQ 1/3 
a= | i 


The separation of the fracture walls is found to be 


gPC 1 -— v2 ) 
ce 


where r is distance from the center of the fracture. The 


(2b e5°7 F242) 


maximum separation B occurs at r=0, and is found to be 


BPC | — vp? ) 
B= a ie Raa tec 8 
1c E 


In order to “use Sores 1..1 fancr egr2.).o) cos cue 





calculation of the radius and maximum separation of the 
fracture formed under the ground surface, it is necessary to 
take account of the effect of the free ground surface. It 
has been, however, found that the influence is only 2-3 


percent if the radius of fracture is less than one half of 
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the depth of fracture, and is almost zero if the radius is 
less than one fifth of the depth (Sun, 1969). 

The fracture model from which Eq.2.1.1 derives ignores 
the effect of cohesive forces near the edge region of a 
fracture on the stress and strain field. The fracture model 
yields an infinite tensile stress at the edge of the 
fracture. A more realistic value for the radius of a 
fracture should be derived from fracture models which take 
account of cohesive forces and finite tensile stress at the 
edge of the fracture. 

One of the modifications of Eq.2.1.1 comes from the 
fracture model proposed by Barenblatt (1962). In his concept 
of a fracture, the fracture is divided into two regions. In 
the inner region (radius a), the opposite faces of the 
fracture are relatively far apart, hence there is _ no 
molecular interaction between them. In the edge region 
(radius ato a'), they are sufficiently close to each other 
so that there are cohesive forces between them. For the 
mathematical simplification of the problem, it is assumed 
that: 1) The width of the edge region of the fracture is 
small when compared with the size of the whole fracture, 
and, 2) When the fracture extends, the shape of the section 
normal to the fracture surface in the edge region does not 
depend on the pressure in the fracture and is always the 
same for a given material under given conditions of 


temperature and overburden pressure. 
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Under the particular case in which the forces of 
cohesion are small enough and may be neglected when compared 
with the overburden pressure and the injection pressure P,, 


the following relation is found. 


OX = a, a (2elee 
Cote a? ee ( We - poh) fete C2 eo) 
where p is density of overburden material, g is 


gravitational acceleration and h is depth. Then, the average 

pressure over the entire circular area of radius a' can be 
2 : 

assumed as Q& P, where P=P,-pgh. In this case, Eq.2.1.1 and 


Eq.2.1.3 are modified in the form 


P ( 3EQ VS 
a’ = oi i ois | So AT ae ve 
hGH Shira? 2) OEP a oa 
8¢ l- yp?) 
B = ———_——— oP Cees) 
CE 


(Sun, 1969, Eq.19 ana Eq.20). 

The theoretical surface deformations of this fracture 
model can be expressed by the product of the maximum 
separation (B of Eq.2.1.3 or Eq.2.1.7) and a quantity which 
depends only on geometric parameters, such as observation 
coordinates, the depth and radius of the fracture. For 
convenience, the fracture model in which the magnitude of 
surface deformations is calculated in terms of B defined by 
Eq.2.1.3 will be called "the first model", and the other 
model related with Eq.2.1.7 will be called "the second 


model". 
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The calculated surface uplifts in Fig.2.1.2 are based 
on -"the second model”. Comparison between the calculated and 
the observed uplifts Supports the validity of the 
mathematical fracture model. Fig.2.1.3b shows that, while 
the cored grout sheet has a nearly circular shape, the 
injected grout moved generally toward the north-east. In the 
other experiment shown in Fig.2.1.3a, the grout sheet has an 
elliptical shape, the movement of the grout sheet is 
greater, and no correlation can be found between ena shape 
of the grout sheet and the pattern of uplift. The maximum 
thickness of the grout sheets measured in cores is much less 
than the calculated values. For example, the theoretical 
value for the case of Fig.2.1.3a is 7.92 cm and the observed 
thicknesses are on the figure. The discrepancy between the 
measured and calculated thicknesses was interpreted as due 
to the fact that the liquid phase of the injected slurry had 
been squeezed out and the soild phase had been compacted by 


the overburden pressure. 


(2) Holzhausen et al. (1980) 

A second example is the initial cycle of steam 
stimulation project of Canada Gulf Resources (1979) 
conducted in the Athabasca oil sands region, Alberta, 
Canada. The hydraulic fracturing was carried out over a 
depth interval of 308-317 m of an injection well in the 
Cretaceous McMurray Formation. The ground deformation was 


measured with a circular array of eight tiltmeters at a 
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radius of 122 m from the well and an additional tiltmeter 
installed 610 m from the well. These nine tiltmeters were 
installed in shallow boreholes less than 7 m deep. The 
tiltmeters measure tilts in two perpendicular vertical 
planes, and have a resolution of 5x10-? microradians and a 
range of 5x10‘ microradians. 

The hydraulic fracturing was initiated three times by 
cold-water injection. The fracture caused by the injection 
was interpreted as vertical because the fracture propagated 
at downhole pressures of about 5.6 MPa whereas the minimum 
possible value of overburden pressure was about 6.0 MPa. 

After the cold-water injection, steam injection took 
place and continued intermittently for about 40 days. Ground 
deformation during this phase occurred episodically over 
periods of a fraction of an hour to a few hours. Gradual, 
long-term changes, resulting from the increasing volume of 
injected fluids and thermal expansion of the formation, were 
reported to prove less distinguishable and less diagonostic 
of formation response than the episodic Neer term changes, 
which were called "events". 

The short-term events were typically preceded by 
gradual increase of wellhead pressure over periods ranging 
from several hours to more than one day. After the wellhead 
pressure rose up to ae peak value (8.0-9.4 MPa), it fell 
rapidly, reaching a minimum (6.7-7.9 MPa) within a few hours 
and then rising again gradually. The tilt response during 


each of the events first changed rapidly as_ pressure 
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dropped, then recovered toward its former value. In most 
cases, the recovery of the tilt field began before that of 
pressure. The overall tilt changes indicated uplift of the 
ground surface during the initial pressure drop, followed by 
subsidence of the surface not in excess of its initial 
uplift. Eight events among eleven detected during the period 
of steam injection were analyzed quantitatively. Figs.2.1.4a 
and b shows the changes of the tilt vectors during the first 
event and during the eighth event, respectively. Figs. 301.5a 
and b show the records of wellhead pressure and the tilt 
changes at an instrument site during the first event and 
during the sixth and the seventh events, respectively. 

The surface deformation during the steam injection 
indicated the occurence of horizontal fractures. This was 
Supported by both the radially symmetric uplift of the 
ground surface and the records of pressure. In any case, the 
recorded pressure did not drop below the estimated 
overburden pressure before the beginning of subsidence of 
the uplifted ground surface. 

The transition from the initial formation of a vertical 
fracture caused by the cold-water injection to the 
subsequent growth of horizontal fractures was accounted for 
in terms of thermal expansion of the oil sands during steam 
injection. "As heating progressed, some vertical expansion 
of the material should have been possible, whereas 
horizontal expansion may have been largely suppressed 


because of the confining effect of the surrounding rock. As 
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temperatures increased, this fracture should have been 
closed by thermal expansion of the oil sands. The resulting 
horizontal aiiperesscon of expansion would lead to an 
increase in the magnitude of the horizontal in-situ stress 
components, leaving the vertical stress largely unaffected. 
It thus seems likely that horizontal fractures began.to form 
when horizontal stresses began to overcome the vertical 
in-situ compression." 

The assumption of impermeable rock does not hold for 
this test site. Therefore, the radius of fracture can not be 
estimated from the relationship between it and the total 
volume of injected fluid, Eq.2.1.1. In fact, it was reported 
that the capacity of the fractures estimated from the 
measured tilts was only 0.2-2.8 percent of the injected 
volume of water-equivalent steam. 

Estimates of the radius of the fractures were based on 
"the first model” for a horizontal penny-shaped fracture. In 
order to estimate the radius from the observed tilt, the 
elastic constants and the driving pressure of the fracture 
should be known. ( This statement holds true for this 
experiment because of the configuration of measurement 
points and the assumption of the shape of fracture.) For 
this experiment Site, the values of the density of 
overburden, the injection pressures at the bottom hole and 
ane elastic constants of the surrounding rocks were 
uncertain. Although the well-head pressures were monitored, 


it has been thought that estimates of driving pressures vary 
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by factors of four or five even when well head pressures are 
measured. 

The uncertainty of the driving pressure of the 
fractures was overcome by uSing the concept of fracture 
intensity factor and fracture toughness. Linear elastic 
fracture mechanics states that a fracture propagates when 
the stress intensity factor at the fracture tip reaches a 
critical value (Kic, called fracture toughness) of the 
material. For a penny-shaped fracture in an infinite, 
isotropic and homogeneous elastic body, the relation between 
the stress intensity factor Ki, the driving pressure Pd and 


the radius a is given by 


K, = ould nee TC (271. 8.) 


G@mccls) OPThe “Seffect = of tthe Mafree-surface onthe, stress 
intensity factor is shown in Section 2.3 ). When the 
fracture is in a state of mechanical equilibrium, the stress 
intensity fractor is postulated to be equal to the fracture 


toughness. Hence the driving pressure is given by 


Gee Au aaa (3.4.9) 


The fracture toughness of the oil sands was assumed to 
be under 330 KPa ./m by taking account of the fracture 
toughness of about 330-1100 KPa,/m for most rock. The 
equilibrium driving pressure just prior to structure 
collapse was deduced as 35-140 KPa. Moreover, for the 


calculation of the radius of the fractures presumed to be 
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horizontal and penny-shaped, the shear modulus of the 
overburden materials was assumed to be between 1.2x10° and 
2.6x10° MPa and the Poisson's ratio to be 0.3. The radius of 
the fractures, respectively distinguished as a short-term 
event, was estimated by using these possible values of the 
elastic constants and the driving pressure, and by comparing 
the magnitude of the measured tilts with the magnitude of 
the theoretical tilts which vary with variation of the 


radius of fracture. 


(3) Evans et al. (1980) 

A third example presents results of a tiltmeter survey 
conducted in 1979 in order to determine the geometric 
characteristics of fractures associated with nitrogen gas 
injection into the Devonian shale in Know Country, Ohio. 
Bedding of all the formations around the injection well is 
approximately horizontal, and the shale contains high angle 
faults of approximate orientation N70°E. 

The injection well penetrated to a depth of 341 m, the 
upper 322 m of which was cased. The casing was perforated 
between 305-314 m. The total volume of the nitrogen injected 
during about 28 minutes is 400 m*® calculated in terms of the 
gas volume at a temperature of 10.2°C and a gas pressure of 
6.205 MPa. Eight tiltmeters were installed in a circular 
array around the well. A shallow watertable was reported to 
have contributed significantly to the background noise level 


of three instruments. 
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The measured tilt changes are shown in Fig.2.1.6a. From 
ene) inspection of the set of tilt waveforms, the tilt event 
was divided into three phases. The tilt vectors accumulated 
during each phase are shown in Figs.2.1.6b,c,d. The 
transition from Phase I to Phase II was identified by a 
sudden change in the previously consistent tilt rates and 
directions recorded by all instruments. The onset of Phase 
II marks a change in the evolution of the fracture growth. 
Wellbore pressure and flow rate remained undisturbed both 
during and following the transition. 

The observed tilt vectors of Phase I suggest the 
occurrence of a long, almost vertical fracture. Based upon a 
two-dimensional model for the free-surface displacements and 
tilts due to an infinitely long fracture of arbitrary dip, 
height (length along the dip direction) and depth to center, 
Subjected to a uniform internal pressure, (Pollard et al., 
1979), the strike, dip-angle and height of the fracture were 
estimated. The complex tilt responses of Phase II were not 
“eee quantitatively. The transition from Phase I to 
Phase II was explained by two possible hypotheses. 

Fig.2.1.6b shows that the tilt vectors are almost 
perpendicular to a direction of about N60°E and distributed 
symmetrically with respect to the strike. The dip angle and 
the length along the dip direction were determined by 
comparing the relative amplitudes of the observed and the 
theoretically predicted tilts. The form of the tilt field is 


dependent only upon relative fracture geometry. Hence, 
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assuming the center of the fracture remains at the depth of 
injection during growth of the fracture, the height and dip 
angle of fracture can be estimated from the relative 
amplitudes of the surface tilts. It was concluded that the 
observed tilt changes of Phase I could be reproduced almost 
exactly by a suitable pressurized fracture of half-length 
183 m having a dip to the northwest 87 degrees. 

The procedure of taking the ratio of the observed tilts 
to the predicted tilts is used to eliminate the uncertainty 
of the other unknown fracture parameters, such as _ driving 
pressure and elastic constants. Assuming the appropriate 
bounds on the shear modulus, 6.9%103= 1702103 MPa... and. the 
Poisson's ratio of 0.25, the bounds on driving pressure were 
estimated as 45-265 KPa. These values are understood to 
represent the possible range of average pressure operating 
in the fracture. Initial formation breakdown occurred under 
modest pumping rates at a downhole pressure of 6.204 MPa, 
which thereafter increased to 8.962 MPa as a result of some 
increase of injection rate. Assuming this breakdown pressure 
is approximately equal to the in-situ stress normal to the 
fracture plane, a simple computation suggests that a driving 
pressure of 2.758 MPa acted in the fracture. It is, however, 
not necessary that such a high pressure extends appreciable 
distances from the wellbore. Indeed, consideration of the 
fracture toughness of rock makes it extremely unlike that 
driving presSures operating near remote fracture edges 


exceed hundreds of KPa. Consequently, the values of driving 
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pressure. deduced from both of the surface tilts and the 
breakdown pressure suggest that pressure gradients exist in 
the fracture. | 

Two hypotheses were presented to explain the transition 
of the observed tilts. First, a comparison with the tilt 
vectors of Phase I and II reveals an approximate reversal in 
trend on instruments located in the northeast quadrant and 
Suggests a partial collapse of the fracture opened during 
the initial phase. This can be anticipated if the fracture 
intersects a natural fracture system at the end of Phase I. 
Also the sharpness of the reversals of the instruments #1, 
#2 and #3 suggests that this intersection would have occured 
to the north-east of the wellbore. However, the tilt vectors 
from instruments #4 and #5 do not conform well to this 
scheme. 

An alternative interpretation is that the fracture 
ceased to grow vertically at the end of Phase I and instead, 
flipped over into the horizontal plane. Physically this 
could occur if there was a major change of the minimum 
principal stress direction, such as that from horizontal to 
vertical. The depth, D, at which the change of the principal 
stress direction is possible was calculated from 
Y.2-4F, 

[—¥Yy 


(22 e.403 





P (2) = P, 00) i. 


where Ph(z) is the horizontal principal stress at depth z, 
zaPvy is the vertical principal stress at depth z, APv is the 


overburden pressure gradient and )Y is Poisson's ratio. At 
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first, Ph(0) is calculated by using the appropriate values 
of p_ , Pv(z)=zaPv (overburden pressure), Ph(z) (breakdown 
pressure) and z (injection depth). Next, the depth D is 
estimated by the relation 

D> ct cle Dele (2.1.11) 

aP, - Py / Cl-¥) 

which is obtained by substitution of Ph(z)=Pv(z) in 
Eq.2.1.10. The calculated depth D for this site is 170 om. 
This value fits well with the result from the analysis of 
Phase I, that is, the center depth of the nearly vertical 
fracture is 335 m, and the half-length of the fracture is 
183 m. It was also pointed out that this estimated depth 


Places the top of the fracture near a lithologic boundary. 


(4) M. D. Wood Inc. (1979) 

A fourth example presents the measurement and analysis 
of surface tilt responses associated with four hydraulic 
fractures produced in horizontally bedded oil shale 
formation and a subsequent injection of explosive slurry 
into these fractures. This field experiment was conducted 
near Rock Spring, Wyoming, in 1978. 

The four hydraulic fractures were produced separately 
at slightly different depths of 124-118 m in a well during 
one month. The explosive slurry was injected through four 
wells, each feeding one fracture, three months after the 
last hydraulic fracturing operation. The period of each of 


the hydraulic fracturing and explosive injection operations 
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was on the order of 20-30 minutes. Eight tiltmeters were 
installed in.a circular array at a radius of 48.8 m from the 
injection well. One instrument was 76.3 m away from the 
injection well. 

The measured tilt changes from the initiation of the 
fluid receion to the shut-in are shown in Fig.2.1.7a for 
me hydraulic fracturing, and in Fig.2.1.7b for the 
explosive injection. The fact that the radial tilt changes, 
in all cases, are non-uniform indicates that the fractures 
are not penny-shaped, but Have irregular boundaries. From 
the sense of the radial tilts, each fracture was regarded as 
horizontal. Whereas the relative magnitudes of the radial 
tilts differs for the four hydraulic fractures, the overall 
forms of the tilt fields are much the same, but are 
considerably different from the form of the tilt field due 
to the explosive injection. 

| To make analysis of the measured tilts tractable, the 
radial oe a changes at eight instrument sites were 
individually assumed to have been produced by a horizontal 
penny-shaped fracture. The result is eight different radii 
for each fracture, which were regarded as an estimate of the 
fracture boundary. The computation of the radius of the 
fractures is based upon the same method as one used in the 
second example cited previously. The value of fracture 
toughness used in the calculations is 373 KPa,/m, and the 
shear modulus and the Poisson's ratio are 3.3x10°* MPa and 


0.16, respectively. 
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Gey Some other examples of ground deformation monitoring 
for the purpose of delineating the form of hydraulic 
fractures are, 
i) Tilt measurements in tunnel for hydraulic fractures 
produced near a tunnel in bedded ash fall tuff at a depth of 
415 m in Nye Country, Nevada, (M.D.Wood et al, 1981) 
ii) Tilt measurements for massive hydraulic fractures at a 
depth of about 1500 m in Wattenburg gas field, Colorado, 
(M.D.Wood et al, 1981). | 
Except for an example cited in (1), the other field 
experiments with tiltmeter surveys were conducted by 
M.D.Wood. Inc. ( The present name is Fracture Technology 


Inc. ), California, on and after 1978. 
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Tables2. 1s lasanGeDiwen. The major content 
of the past field experiments in which 
surface deformation measurements were used 
to investigate hydraulic fractures. 
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2.2 Inception of Hydraulic Fracture 

The inception of hydraulic fracture has been mainly 
examined from the stand point of the occurence of tensile 
failure on the wall of wellbore. There is also a possibility 
that the hydraulic fracture breaks out as shear failure. The 
wall of wellbore or the surrounding formation must contain 
some pre-existing cracks. Some wells must allow the injected 
fluid to flow into the formation. Circumstances of 
well-fracturing can not be limited to a stereotyped form. 
Various criteria have been proposed regarding fracture 
initiation from each of these stand points. Here I will 


review some of those criteria. 


(I) Criterion based on tensile strength of smoothed well 

The criterion proposed first by Hubbert et al.(1957) is 
that fracturing will start at points on the wall of wellbore 
where the Regncri ve stress is equal to or greater than the 
tensile strength of the rock. Based on this criterion, we 
can predict the direction ‘and critical pressure of the 
fracture inception, if the circumstances around the wellbore 
are consistent with simple models. In order to understand 
the stress field around an imaectiba well, we have to take 
into consideration at least four kinds of effects. These 
effects arise from (1) the regional tectonic stresses, (2) 
the injection pressure, (3) the formation fluid pressure and 
(4) the fluid flow through the porous rock from the wellbore 


into the formation. In the following, it is assumed that 
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rock around the wellbore is infinite, elastic, homogeneous 
and isotropic, and is under a nonhydrostatic state of 
regional stresses with one of the regional principal 
stresses acting parallel to the vertical axis of the 
wellbore. 

First, the wellbore is approximated by an infinitely 
long vertical cylinder of radius a, in a very large 
rectangular plate of infinite thickness on which edges the 
horizontal stresses 6, and GO, act, as shown in Fig.2.2.1. 

The stress field around the cylindrical hole due to the 


horizontal compressive stresses, Ox, and Oy , is given by 





2 2 : 4 
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where compression is considered as positive. 
Next, let us consider the stress field around a 
cylindrical ring subjected to the internal pressure P 


(Fig.2.2.2). The stress field with the boundary conditions 








Ss CheekeeP) 9 8, Cr baeeo (2.223) 
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PAO. A a ie ® Gy. : A cylindrical cavity 
subjected to far-field stresses. 





ee Pe eae ae A cylindrical cavity 
subjected to an internal pressure. 
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Ser = 0 
The total stress field around the hole subjected to the 
internal pressure P and the far-field horizontal principal 
stresses OG, and G) is obtained by superposing Eqs.2.2.1 and 
2.2.3. When the outer radius b is considered to be very 
large in comparison with the inner’ radius a, the total 
stresses at the wall of wellbore, r=a, are expressed in the 


Simple form 
qr) he ary (2) 


is Sij + $ij 
Tre = P (oro ea5 
hes = (6+ G )-2(6%- 6 ) cos 26- P 
ibe Eee. 


Here we assume that fracturing will first occur at 
points on the boundary of wellbore where some stress is 
equal to or greater than the tensile strength To of the 
rock. Eq.2.2.4 indicates that the only stress which can 
become tensile at r=a is the tangential component, Tos, and 
that the point which stress will first become tensile is 
cos26=-1.(62%+T/2) , and the stress is given by 

3 6x - inten G 

where =6, <,G, has been assumed... The, fracture will start 
along a plane normal to the direction of the minimum 
Principal stress. Denoting the critical pressure which 
overcomes the tensile strength of the rock by Pc, we have a 
Criterion of the fracture initiation 

Fogo ee toe ecs! 
(Hubbert et al., 1957). 
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Eq.2.2.5 is revised when the presence of pore pressure 
is taken into consideration. The difference between the 
internal pressure P and the initial pore pressure P, will be 
denoted by 

ei aed ed a OL202 66) 
and the effective stresses of 6, and Gy by 

Oxx : Ox 5 A 

Syy = OF - F, 
(cf., in general, Eq.2.2.7 is given by Go 


(2727) 
Sete jf ~k SP, 7 


Then the effective stresses around the cea in porous media 
Of pore pressure P,, caused by the far-field stresses 6G, and 
OF and the internal pressure P, are given by 

(2e) cle) 


ai kor ie t. Sij 


where Sij (“© is obtained by replacing o, and G, in Eq.2.2.1 


(2e) 


(2.2.8) 


with as is Gyy, and Sij*’is obtained by replacing P in 
Bones exo with Pw. Consequently, the criterion of the 
fracture initiation is given by 
(2e2i¢2) 
or meni eae Oe ee. 
(Haimson et al., 1967,  Eq.12), where Pe is the critical 
value of Pw. As shown tacitly in the above discussions, as 
long as such an infinitely long vertical cylindrical hole is 
used as a model of well-fracturing, there is no prediction 
of the occurrence of horizontal fracture. 
In an open-hole well-fracturing procedure, the 
prospective producing horizon is packed off and fluid is 


introduced into this zone. Kehle (1964) modelled the 
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circumstances of well-fracturing by a band of uniform 
pressure and two bands of uniform shear stress acting ina 
cylindrical cavity as shown in Fig.2.2.3, and obtained the 
conclusion, from numerical analyses in which the radius of 
the hole was taken to be 0.21 m; that, for any packer length 
(1.5, 3.0 and 4.6 m) and any pressure interval (from 1.5 to 
15.2 m at 1.52 m increments), the magnitude of the induced 
tangential and vertical stresses due to the shear bands is 
negligible (at most a few percent) in comparison with those 
caused by the band of pressure. Fig.2.2.4 shows the result 
for the case of packer length 3.0m. 

Two interesting regions of induced stress are, (1) 
either end of the pressurized interval where the tangential 
stress is zero and the vertical stress is approximately 94 
percent of the pressure, and (2) the central part ( about 80 
percent ) of the packed-off interval where the tangential 
stress equals the pressure and the vertical stress is aoe 
The induced stresses in the central region of the packed-off 
interval is shown to be almost identical to those predicted 
when the cylindrical hole is loaded along its entire length. 

| The induced stresses on the wall of wellbore, 
are summarized as follows. 


In the central region of the packed-off interval, 


Be at in San aes! ot ER Son = 0 a LEY 
in the very ends of the interval 
So) -tipe, tec eteepses , fins 21.208 94. P (2.2.10b) 
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Pigede2. 3 a... Schematic diagram of ( 
hydraulic fracturing set-up and ( 


associated boundary stress assumed y 
Kehle. 
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Fig.2.2.4....«., Ventical..and..tangential 
stresses induced at the wall of the hole 
dia tomes tresses) in. Fige2ee.>c. | - from 
Kehle, 1964 ). 
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The total Stresses at r=a, Tij , are given by 


Superposing Eq.2.2.1 and Eq.2.2.10. 
ate eek. Sh ¢3> C2eoinl | 
Ff 4 = 555 + Siy ) 
When the pore pressure P, iS considered, the effective 


stresses of Tj are expressed in the form 
c3e) cle) mse) 
5 x Si; a Si5 (2 eae) 


where aie, is obtained by replacing P in Eq.2.2.10 with Pw 


of Eq.2.2.6. The concrete expression of hie a se 


in the central region, 


(3e) 
Tet oud 


c3e) 


Doet Fase at Theo SyuTa (s‘acbore Se +250) 


SSH — OR SS Poy -*2E, 


J 
Tee on ite e - P, (202. 13a } 
and in the end region, 
1 ee P, 
Giect? 3 Ok Ovy (ae eo stn) 2.) 
syn = Sy ca 2k, 
Tos = Gyy — 0-94 Py = Gy- Py - 0-940 P-P,) (2.2.13b) 


(Kehle, Eq.37), where 63, is the effective stress of the 
vertical component SG, of the regional stresses. 

From Eq.2.2.13, we observe that, besides the tangential 
stress Jp, in the central region, the vertical stress T,,°” 
in the edge region can become tensile. This suggests a 
possibility that horizontal fracturing takes place. The 


critical pressure Pc for the occurrence of a horizontal 


fracture is given by 
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For 422 or > — ~ - (2.2 1840) 
0.94 0.94 


Next we will suppose that fluid flows from the borehole 
into the surrounding formation. To simplify the flow 
problem, the wellbore is again assumed to be an infinitely 
long cylindrical cavity.: It is also assumed that the 
formation has uniform permeability so that the fluid flow is 
axisymmetric, and that the injection fluid has properties 
Similar to those of the formation pore fluid, and that the 
fluid flow obeys Darcy's law. 

The effect of the fluid penetration is demonstrated by 


the stresses rhe 


caused by the pressure distribution P(r) 
in the formation under the condition of steady flow and with 
the boundary conditions 
ca) dl” Prk SS ale a 
Berea) = s5 (tb), =70 Sear TEE §, 
med) PP (rsb) =. 0 eh Ae ng fot ty 


The stress tensor Sie due to the pressure P(r) is given by, 


in terms of the radial displacement u, - 


(4) du 
sven” A4+ <j eres (3+ 2 xP 
C4) Uu 
= trot ows ( SX 2) Ot 
Summa. eM cta2 Se pap 
43 (2.2600 
5 saisi a 4 cores) P 
du u 
dr r 


and other stress components are zero, where \ and Kare 


Lame's constants, and & is a material constant. Since we 
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are supposing the state of steady flow, the pressure 
distribution P(r) satisfies the Laplace equation. The 
boundary . conditions Eq.2.2.15b determine a solution of this 
differential equation. 

PCr) = P log(b/a) / log(b/r) ‘psy ig a 
On the other hand, the boundary conditions Eq.2.2.15a 
determine a solution of a differential equation in terms of 
u, which derives from the equilibrium equation. Finally, the 


stresses Si are given by 








ri 
(4) XEP b a’ b B 
Spr 2(I-v) log (b/a) 3 r b?- a2 d a 
z v. 
OS aEP jee eae: tog} 
as 2(1-Y¥)l03(b/a) Rereibe= a. E z 
2 b 
Secor aEP {1- pfoqee og — | (222082 
2(i-V) 109 (b/a) fe vasa a 


(cf., Takeuchi, P149). When b is considered to be very large 
compared with a, the stresses Eq.2.2.18, at r=a, are 


Simplified in the form 


aii: 0 
<4) ae a 
2 careers 
See fa) oaearsg 
Peed 
(4) _ 
Yes KE Breas FP 


The total stress field is obtained by the superposition 
s 
of the stresses Sue due to the far-field regional stresses 
Gx and Gy , Ss et due to the pressure difference Pw, 


Si et due to the fluid penetration and the initial pore 
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pressure P,. 


ae) 


(4? 
t: ~ Sij Si; + Si j +P (2.2520) 
For convenience, let 
e i= 2) ( y) > Cy 
OE ———""= ACi-2Y) = = = ya 
bees B Wee 75 Cy ee eee 


(cf£., Rice et al., 1976, p239) where Cr and Cb are the rock 
Matrix compressibility and the rock bulk compressibility, 


respectively. Then the stresses aire. are given by 


ee sD ae 
C4) _ . e 
Teg = -§ Oxx- Sy yJ-2( 6,4 S,, )cos 26-P, + 2nP + P. 
4) 
(a eR (Ae esas (7-3) 23} 
From Eq.2.2.22, the minimum value of T,,*? is found at 
8 = Le 1/2, 
@,min _ sf - 
1h = Sips Oo ac tld Coon cea) 
The effective stress of cen ee is 
(,min 
fan. a Creer = 5 OX. - Syy- 2A +2NP, (2.2424) 


The critical value Pc of the pressure difference Pw, which 


will break the wall of wellbore, is now given by 


Py ate (2.2.25) 
2(!-n) 
(Haimson et al., 1967, Eq.10). 

When the induced vertical stress T,,* in Eq.2.2.19 is 
combined with T,,°° of the Kehle's model (Eq.2.2.13b), the 
critical pressure of the occurence of horizontal fracture is 
given by 
To* Sz 
cones 1.94-21 


Vv 


(202526) 
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As o<A<1 where @ is the porosity of rock, and 
Opcsel erese for-e0rock, *1t *f6éllows ONE K Pf eFSThe 
comparisons between Eqs.2.2.25 and 2.2.9 for the horizontal 
fracture and between Eqs.2.2.26 and 2.2.14 for the vertical 
fracture suggest that the fluid penetration causes the 
inception of either type of the fractures at lower injection 
pressure than that for the non-penetration case. It is. also 
Suggested that very rapid pressurization leads to a higher 
fracture pressure than needed when pressurized fluid is 
allowed sufficient time to penetrate the walls of the 
wellbore. 

Heat conduction through solids gives rise to stresses 
and displacements analogous to those caused by fluid flow 
through porous media. Therefore, the above discussions can 
be applied to estimate thermoelastic effects associated 
with, for example, hot fluid injections. The stress field 
due to non-stationary fluid flow has been studied by Rice et 


al.(1976) in more general form. 


(II) Criterion based on fracture toughness of notched well 
The previous criterion is based on a critical tensile 
Stress fracture and neglect of consideration of cracks in 
the wall of the wellbore. Here, I will refer to a criterion 
of fracture initiation which is based on linear elastic 
fracture mechanics and allows pre-existing cracks on the 


wellbore. 
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Assume a cylindrical hole of radius a with two radially 
opposed cracks of length L subjected to an internal pressure 
Pand the far-field stresses G7 and 6, , as .shown in 
Fig.2.2.5. The stress intensity factor Ki at the crack tip 
is given by 

K, = PYRL F(L/a)~(o; cos*x +65 sin% )F(L/a Wi 
+ (G) cos2k-6,¢0S82%) G(L/a) VTL (2227) 
(Abou-Sayed et al., 1978, Eq.10), where F(L/a) and G(L/a) 
are shown in Table 2.2.1. Let Kic and Pc respectively denote 
the fracture toughness and the critical pressure for crack 
advance. From Eq.2.2.26, the critical condition of the crack 


extension is expressed in the form 


Rie 
P. = 1(&) (67-63) +03+ ——_—_— 
F(L/a) wo 
. G(L/a2) 
I(a) = ¢cos*X— —— cos2KX Creo) 
F(L/a) 


The function I(*) is shown in Fig.2.2.6. Eq.2.2.28 states 
several features of the fracture initiation under’ the 
circumstances of the wellbore with such pre-existing small 
cracks. 

bi Wiy On , the orientation of the pre-existing 
crack does not contribute to the critical pressure, and only 
the length of crack influences the critical pressure by a 
factor of Keer FcL/a)dn , suggesting that fracturing in 
the Peer with larger cracks will start at lower 


pressure. If I(~)= 0 , the critical pressure depends on 


either principal stress. 
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a's ee Pts hea ae On Crack in a borehole wall 
subjected to far-field stresses and 
internal pressure. 
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E100 Se Galen s. 1(X) defined by Eq.2.2.28 
for different ratio of L/a. ( from 


Abou-Sayed et al., 1978 ). 
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If x = 0 F Vie u, the orientation of the 
pre-existing cracks is perpendicular to the direction of the 
minimum principal stress axis, the crack advances most 
easily. If & = H/2 , the maximum pressure is required to 
advance the erackjcibecause I[fd)p a I(AdiniXo2 B from 


Fig.2.2.6. When the case of L/a=0.2 is examplified, the 


difference of the critical pressures for & =0 and 
Md=OR/2 is given by eS 46 Sytt0Gs ) ; where 
L(cd 30/2 )4-4.25 and I(&=0)=2=-0.25 are read from 


Bragee.270.200 thesrother.-hand, .if;«the.~.crack length..is 
considered to be large enough compared with the radius of 
wellbore,@a,i.e., L/a ~ @& , ne follows that 
1(m = 0) ~ [(m&=7/2), hence there is no difference between 
these critical pressures. 
BE = Om andes i= = Oi o/isay, L/a=0.1), 

2P. ~ 365-9, +2K,. /2. 264m (242629) 
where the value 2.26 comes from Table 2.2.1. It is observed 
Phat inbCwsalNesEQi2a2-29) isjrelated.to..363-6, by.factor,2, 
whereas in Bq.2.2.9 the factor is 1, and that &q.2.2.29 
holds true without regard to the existence of pore pressure. 
It is noted that Eqs.2.2.27-29, in principle, can not used 
for the case L=0 and there is no consideration of the stress 


intensity factor for shearing mode. 


(III) Mohr-Coulomb criterion for shear failure 
Coulomb (1733) proposed a criterion that when. shear 


failure takes place across a plane, the normal stress and 
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the shear stress across this plane are related by 

oh iS So + wa eee eu 
where S, is a constant which is regarded as the inherent 
shear strength of the material, and HW; is a constant which 
is called the coefficient of internal friction of the 
material. The coefficient W; is usually expressed as 

y. = tang (202650) 
and @ is called the angle of internal friction. 

The two-dimensional case will be considered. The normal 
Stress © and shear stress CT acting across a plane 
perpendicular to the 0©,~0, -plane and making an arbitrary 
angle & with the direction of least principal stress 0, , 


(Fig.2.2.7a), are 


26 = (0,+65) + (6,-635)cos2m 
: (2.2.92) 
2T = (6,-G3 )Sstn2X 
A very convenient method of graphically representing these 
expressions, known as the Mohr diagram, is obtained by 
plotting values of normal and shear stress with respect to 
G@ - © -coordinate axes for all possible values of the angle 
eee. (Pagc2.2.7b). “Then the- relation’ of Eqs2.2/30" “is 
represented by a Straight line, for example, AB in 
ie ape ieee aie B= 
Eq.2.2.30 suggests that failure will not take place if 
the values of O and T so found lie below the line AB, and 
failure) (will) (take “piace if the ‘circle on dc; and Gy as 


diameter just touches AB. The substitution of Eq.2.2.32 into 


Eq 202730 gives another expression for this failure 
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criterion 
[Cw 2+ Fp] = 65 [Cu,7 +1)? Hy] = 28, 
{ 


Pe e33) 
at & = +(T/4to) 


7é., cot2X = ft 2ui 
When the presence of pore pressure P, is taken into 
consideration, Eq.2.2.30 is revised as 
[tl 2 So+H; (O-P,) (2.2.34) 
As suggested by a shift of the Mohr circle in the left 
direction along O -axis in the Mohr diagram, an increase of 
the pore pressure P, accelerates the occurence of failure. 
Mohr (1900) replaced the criterion E£q.2.2.30 by a 
functional relation characteristic of material 
[tl = f(r) (PIAS) 
The curve of (0) is called the Mohr envelope, and obtained 
experimentally as the envelope of the Mohr circles 
corresponding to failure under a variety of conditions. With 
regard to the function f(0) , the coefficient of internal 


friction is defined by 


= ale paben pe i ees G 


Tc 

ies: $a 

The angle of internal friction has values usually between 

20-50 degrees and most commonly not far from 30 degudss ( 

Bi ~ 0.6). The Mohr envelope is usually concave downward. 

At sufficiently high pressures, nearly all rocks deform 

plastically, and the Mohr envelopes becomes approximately 
parallel to the O -axis. 

The tangential and normal stresses acting on a plane 


with an arbitrary orientation relative to three principal 
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Stress axeS can easily be seen from the Mohr's 
representation of stress in three dimensions (Fig.2.2.8). As 
long as the failure strength through the body of the rock of 
interest is assumed to be characterized by unique values of 
So. and i in Eq.2.2.30 or a unique function f in Eq.2.2.35, 
the Mohr-Coulomb criterion basically states that (1) 
possible failure planes are always parallel to the axis of 
the intermediate principal stress, and (2) this stress does 
not affect the angle of the failure plane. It is so because 
epee line |tj=S5,+.9) org the. curve |It)= (0) will first 
touch a point on the circle on 6, and O; as diameter, and 
this circle represents a plane of P=0. Past laboratory 
experiments (e.g., Mogi, 1971; Handin et al., 1967), 
however, have pointed out that the above statement (2) does 
not necessarily hold true. This will be described later. 
Most sedimentary and metamorphic rock is anisotropic. 
Representing the anisotropy by the presence of a plane ee 
weakness, its effect on strength can be seen by using the 
Mohr's diagram. The inherent shear strength and the 
coefficient of internal friction in the place of weakness 
will be denoted by Sw and My, respectively. In Fig.2.2.9, 
let the point P, the line AB t = SwtP,®@ and the line CD 
ton og ees respectively be a stress state on the plane 
of weakness under the principal stresses GO, , G6, and O;, 
the critical condition of the. failure in the plane of 
weakness and that through the body of the material. Suppose, 
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now, that the principal stresses are changed up to o , G 
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The Mohr diagram in 3-D. 
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and OG; as shown in Fig.2.2.9a. The point P will move to the 
point P', and will touch the line AB. Failure first becomes 
possible in the plane of weakness. On the other hand, there 
is another possibility that failure first occurs through the 
body regardless of the presence of the weakness, as a result 
of, for example, the changes of three principal stresses as 
shown in Fig.2.2.9b. 

Extremity of anisotropy is represented by the presence 
of a discontinuity in the material such as a joint, crack or 
fault. Byerlee (1977) points out that the Erictaons | 


strength can be expressed by the relation 


n 


aw 
i 


0.850 : 0-5 <& < 2 kbar 
0.5+0:60 «+: 2.0 < 6 < 20 kbar 


Crease 


and that this relation almost holds without regard to types 
of Bock, temperature, presence of gouge between 
discontinuity surfaces, and roughness of the surfaces, i.e., 
regardless of smoothed surfaces, pre-made shear failure 
Surfaces or interlocking surfaces made by tensile forces. 
The frictional strength under the presence of pore pressure 

P, is expressed by the relation 
T = 0.1+0.6(5=P,) (2.2.38) 

(Byerlee, 1967). 

Criteria of shear failure may be expressed in the 
general form 
‘ F(S, 503,65) = 0 (2.2.39) 
One of the failure criteria which take account of the effect 


of the intermediate principal stress is the von Mises 
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criterion proposed as a yield criterion, and is expressed by 
(6, -6,)°+(5,- 6 )*+ (65-6)? = 20, (> 240) 
where GO, is a value characteristic of material. The von 
Mises criterion has the physical interpretation that yield 
occurs when the strain energy reaches a value characteristic 
of material. On the other hand, it has been shown that the 
relation 
4} (6,- 65)" + (6-55) + (6-7 F = B0+55) (2.2.41) 
explains rock failure well, and the function g is regarded 
to be a monotonically increasing and almost linear function, 
(Mogi, 1971). From the facts that the planes of failure are 
observed always parallel to the intermediate principal 
Stress axis and that (0 ,+03)/2 is related to the mean 
stress acting on a failure plane, this failure criterion is 
interpreted such that failure will occur when the shear 
strain energy reaches a critical value which varies with the 
mean Stress. Past observations (e.g., Mogi,1972) also 
indicate that the angle of the failure plane relative to G@ 
or O; axis depends on the magnitude of O, , whereas the 
Mohr-Coulomb ioritterzon can Mmoti-predict: ther effect of °O,e¢ on 


the failure. 


p 


se Beeoqoiq not 


ee 


2135 dike dus? te 













ag? 


aa 


5 el 3 sisiw 
fi 


nofaedias” ese Ae 


LB 


i> nedw 87B93¢ 


,lelistema 4 


-. 


norse [94 
. & 
m “sy 
\ 6 ae 
‘\F 


2507 
C2Oon 


2ixe i ees 


prijos eaeme 


2 heteiqss 


yotene nse 





2.3 Fracture Orientation Changes 

Two examples of the change of fracture orientation have 
been seen in Section 2.1-(2) and (3). It might be a rather 
common phenomenon that a fracture changes its orientation in 
the process of growth. Some factors of fracture orientation 
changes are reviewed here. 

Holzhausen et al. (1980) reported the initial formation 
of a vertical fracture caused by cold-water injection and 
the subsequent growth of horizontal fractures during steam 
injections. The orientation charge was interpreted as a 
result of thermal expansion of the surrounding materials 
during the steam injections, which led to the build-up of 
horizontal stresses around the initial vertical fracture and 
resulted in a stress state where the horizontal stresses 
exceeded the vertical stress compression. 

Without the thermal effect, change of fracture 
orientation can be explained from the more general 
Standpoint of a change in the minimum principal stress 
direction around an initial fracture due to subsequent fluid 
injections. Dusseault (1980) explained the change .of 
fracture orientation as follows. "IE the injection rate is 
relatively rapid, then energy dissipation due to viscous 
traction on the walls of the fracture can result in thick 
fracture generation with rapid local overstressing, 
particularly near the injection point. It is believed that 
the minimum principal stress is increased to a value 


somewhat higher than the intermediate principal stress 
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before reinjection takes place. Fractures thereafter tend to 
recur by fluid injection in the direction approximately 
orthogonal to the original plane." The Amecress analyses of 
Dusseault et al. (1980) indicate more clearly that regional 
stress changes can be brought about by massive injection 
volume, change of fracture orientation is extremely likely, 
and thermal input enhances this probability. The numerical 
analyses are based on a vertical fracture model in layered 
media with the assumptions of reasonable material property 
and geometry related to the application of hydraulic 
fracturing in the Cold Lake oil sands region. 

Evans et al. (1980) reported that a significant change 
of the surface tilt responses during the continuous 
injection operation suggested a possibility of the 
transition of the fracture orientation from vertical to 
horizontal. The transition was interpreted in connection 
with the in-situ minimum principal axis direction varying 
with the depth. The top depth of the initial near-vertical 
fracture calculated from the measured tilt responses 
coincided with a depth estimated by the analysis of the 
in-situ overburden ee gradient, i.e., the depth where 
the minimum principal stress direction changes from 
horizontal to vertical. Although it can not be simply 
concluded that the latter depth equals just the depth of the 
occurrence of the horizontal fracture, as a rule the spatial 
variation of pre-existing stress field could be a major 


factor in fracture orientation changes. 
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Lithologic factors also must control the propagation of 
fracture. Evans et al. (1980) remarked that the calculated 
top depth of the vertical fracture corresponds to a 
lithologic boundary, and the overlying lithology must have 
acted a barrier to prevent the propagation of the vertical 
fracture. Dusseault (1980) described as follows. "It is said 
for the Athabasca oil sands region that, in those areas 
where horizontal fractures are originally created, the 
fractures tend to climb upward at angles of 10 to perhaps as 
high as 25 degrees. The tendency for a horizontal fracture 
to climb has often been explained in terms of density; if 
air or gas is injected, it is assumed that there is a 
tendency for fractures to climb upward because of the 
buoyancy of the air or gas. Another explanation exists: ina 
material with variably oriented anisotropies (bedding 
planes, joints), fractures will tend to climb along 
appropriately oriented bedding features and other lithologic 


discontinuities.” Even if the uncertain lithologic factors 
are not considered, the tendency for a horizontal fracture 
to turn upward can be predicted theoretically, as shown 
below. 

Several workers in fracture mechanics have proposed 
theories which incorporate the magnitudes of stress 
intensity factors at crack tip, Ki (for open mode) and Kii 
(for shearing mode), and predict the propagation direction 


(e.g., Sih, 1972, Erdorgan and Aksogan, 1974). These 


theories state that a non-zero value of Kii should result in 
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propagation out of the original plane of the fracture in 
directions governed by the sign of Kii. In the following, I 
will refer to two of these theoretical predictions. One of 
them observes the non-zero Kii resulting from the mechanical 
interaction between the free-surface anda fracture. The 
other analyzes the stress intensity factors of a fracture 
Subjected to the internal pressure and non-hydrostatic 
regional stresses. 

Pollard et al. (1979) used a two-dimensional fracture 
model, and analyzed the stress intensity factors at the 
crack tips of an internally pressurized plane fracture under 
hydrostatic far-field stresses in a semi-infinite region. 
The analyses led to the prediction that a shallow fracture, 
whose original plane is not vertical, tends to turn upward 
at the ends and propagates toward the free-surface. 

Assume an infinitely long fracture in an infinite 
region. The stress intensity factors at the tips of this 
fracture are given by 

penne aay sis) 4 Ki 8= 0 CE ESe. 
where P is the uniform internal pressure anda is the 
half-length of the fracture. If the depth of the fracture, 
d, is sufficiently large compared with the half-length, the 
fracture may be assumed to be located in an infinite region. 
Such a fracture should propagate in its own plane, because 
Kii=0. 
For shallow fractures, the mechanical interaction with 


the free-surface must not be negligible. Fig.2.3.1 shows the 
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stress intensity factors normalized by Ke. and plotted versus 
the ratio d/a. In the case of vertical fractures 
(Fig.2.3.1a), symmetry about the long axis of the fracture 
makes Kii zero, but the values of Ki increase sharply as d/a 
decreases to a limit of 1.0. Because of Kii=0, the vertical 
fracture should propagate in its own vertical plane. The 
fact that the value of Ki at the upper end of the vertical 
fracture exceeds that at the lower end suggests that the 
fracture preferentially propagates upward. In the case of 
horizontal fractures (Fig.2.3.1b), the values of Ki at both 
ends are equal and sharply increase as d/a decreases toward 
zero. Non-zero Kii are now observed. The values of Kii at 
both ends are equal but opposite in sign, and sharply 
increase as d/a decreases. The existence of non-zero Kii and 
the signs of Kii at both ends suggest that a shallow 
horizontal fracture should turn upward at both ends and 
propagate toward the free-surface (Fig.2.3.2). It has been 
also observed that varying the inclination, at constant d/a, 
does not result in large variation in stress intensity 
factors, except where the upper end of fracture is very near 
the free-surface. An inclined fracture also will proceed out 
of the original plane, generally tending toward the 
free-surface at both ends, as suggested by the non-zero Kii 
and the signs. 

The comparison of the magnitudes of Ki in Fig.2.3.1a 
and b indicates that, for comparable d/a, Ki for a 


horizontal fracture is greater than that for a vertical 
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Rtg 43.34 209 . 9F.. Interpretation of signs of 
stress intensity factors and propagation 
directions of a fracture. 
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fracture, except where the depth is. very close to the 
Surface. 
| An assumption of uniform pressure acting on fracture 
walls is more or less at variance with reality, especially 
when the frictional interaction between viscous fluid and 
fracture walls is considered (e.g., McClain, 1969). 
Fig.2.3.3 shows the difference of the stress intensity 
factors resulting from a type of pressure gradient ina 
fracture from those resulting from uniform pressurization. 
This figure indicates that the magnitudes of Ki and Kii for 
a fracture characterized by the symmetric linear pressure 
gradient are less than those for a fracture characterized by 
the uniform pressure. On the other hand, it has_ been 
observed that the ratio Kii/Ki for the former model is 
greater than that for the latter model. This suggests that a 
horizontal fracture maintained by the pressure with such a 
gradient should tend to propagate toward the ground at 
higher angles than those in which the pressure is 
essentially uniform. 
When a symmetric, linear distribution of driving 
pressure is given by 
P(x). Py 2 AP [els Pla) =:0 5.0 slxiea (203.2) 
the stress intensity factors for a deeply buried fracture 
(i.e., the fracture may be assumed to be in an infinite 
region) are given by 
WaieeGirs —i2asPiniyaia 40 (KX ,,)-=10 eree 


(Lachenbruch, 1961, Pollard, 1976). In the case in which P, 
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uniform pressure (solid lines) and linear 


(broken lines) distributions. 


from Pollard et al., 
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poeecaken tO oe | twice P in Eq.2.3.1, the total force 
Operating in the fracture is same as that for the model of 
uniform pressurization. Simple computation indicates that 
Bie ratio BO.2.4,3/20.2.3.1° in terms of Ki is about 0.726. 

The stress intensity factor at the tips of a horizontal 
penny-shaped fracture is shown in Table 2.3.1, .in which Ko 
is the stress intensity factor ( Mode I ) for a penny-shaped 
Seecture if am infanite regryon,( cf.,vSection 2.1 ). 

Next, it is shown that fluid injection into a wellbore 
with a pre-existing crack of an arbitrary orientation 
relative to the direction of far-field stresses tends to 
make the cracks extend in a direction which is more nearly 
perpendicular to the direction of the minimum principal 
stress than the original plane was. This is predicted by 
Strain energy release rate hypothesis (cf., Abou-Sayed et 
al., 1978). Consider a crack of length 2L oriented at an 
arbitrary angle & relative to the direction of the minimum 
principal stress, as shown in Fig.2.3.4. The stress 
intensity factors Ki and Kii for the pre-existing crack 
Subjected to the internal pressure P are given by 

K, = (WL *(P-G& sin'x- 63 c0s’X ) (2.3.4) 
Ks CTL '/2+( 6-64) sin 2K (2.3.5) 


(Rice, 1968), where G6, and G, are the maximum and minimum 


compressive far-field principal stresses, respectively. If 
e 


the crack extends in an arbitrary direction Y , the 


corresponding release rate |(¥)is given by 





1.2042 0.0899 
1.1863 0.0822 
Pe tetee 0.0640 
1.1320 0.0514 
1.0759 0.0248 
1.0275 0.0082 


1.00399 0.00062 





a=-radius of fracture 
d=depth of fracture 


Ke=+PVa » P=internal pressure 
iy eB see BS 1 A ae Stress intensity factors 


for a horizontal penny-shaped fracture 
(from Kassir and Sih, 1975, p225). 
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EY FB OS SV ee Orientation of an 
extending crack which maximizes energy 
release rate (from Abou- Sayed et 
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eee et al., 1974), where E and p denote Young's 
modulus and Poisson's ratio of the material, respectively. 
The crack will be assumed to advance in such a direction 
that I'(1) becomes a maximum. Then the relationship between 
the direction of crack advance Tmax and the ratio K,,/K, is 
given by Fig.2.3.5. 

SUS UCLOUMeGteEOG~t.s<F And 2.355 Int” 7H. 2.5.6 
reveals that |’ is proportional to the crack length. 
Therefore, if the crack advances at a prescribed finite 
value of | Preset the conditions for such a_ sufficiently 


long open crack to be stationary are satisfied only when 


)’cr)=o fevice f= tr ~ 2 Aree ae ty 
From Eqs.2.3.4 and 2.3.5, this condition is expressed as 

P = CO sind + Oy cosh (23a 
and OS eros sim 20 = O (2-33.19) 
If O, and OO; are not equal, Eq.2.3.9 implies that 

a = 0 Ore ost TC / 2 
The former case, X& =O  , corresponds to the situation 


that the pressure P equals the minimum compressive stress 
G3 , whereas the latter case , & =+M/2, corresponds to 
the situation that the pressure p is equal to the maximum 


compressive stress O,. 


. 


cu 


Hf 


" 


ti» 


(£9 


st) 7 “ ) 


jon 918 Z0\, brs 


ye] 

> 

» “y 
50 4 
- a 
4 -* 





. 





Ss 


80 


Now consider a problem; to which direction of these two 
Stationary conditions does the pre-existing crack of neither 
&=O nor KX=t1/2 eetand more likely ? We are considering a 
crack tending to open due to the internal pressure, hence 
eee O 8. Since 6, 7 Ox , Eq.2.3.5 indicates that Kii is 
positive for & in the interval ( o, W/2), whereas Kii is 
negative for & in the interval (-%/2,0). This observation 
and Fig.2.3.5 indicate that for & in the interval (0, T/2), 
the predicted angle YT of the extending crack is positive, 
whereas Y is negative for & in the interval (-"/2, 0). 
Consequently, this maximum energy release rate hypothesis 
predicts that the pre-existing crack tends to extend ina 
direction which is more nearly perpendicular to the 


direction of the minimum compressive stress. 
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3. Theoretical Size of Surface Deformations 
In the following, a fracture is assumed to be produced in an 


isotropic, homogeneous, Semi-infinite elastic medium. 


3.1 A horizontal penny-shaped fracture 

Analytical expressions for the surface deformation due 
to a horizontal penny-shaped fracture have been obtained by 
Ren Jen Sun (1969). In this fracture model, it is assumed 
that the surface of fracture is maintained by uniform 
pressure of injected fluid, and that fluid flow into the 
rock adjacent to the fracture is so small that the rock can 
be regarded as impermeable. 

mhe surface deformations at distance r from the center 
of a fracture at depth d and of radius a in a semi-infinite 


elastic medium of Young's modulus E, Poisson's ratio p are 


given by 
Uplift = 
U = f-a- ({ksin$-sScos$) (Seb. 14) 


Horizontal displacement = 


; atalk si ne 
R = T° da: ES SS ee a a a re 
(dtavk cos$) + (atark si n$)? 





dikcos$-adk $ inStrakeos$ 
(dik coss-adksin$+akcos6)* + (ak cos$+diksing tak sin®)” 


(30a. 2) 
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Tilt = 
_ du 
T= ac ec eba Bri 
Horizontal strain = 
dR 
s = ar (aa) 


where 


32 32 (321745) 
2 Aes 2Zad 
seta" ( ratte at ) (3.1.6) 
8(1—-v2) 
rescuer onal (3.1.7) 


and P is the internal pressure of fracture (Ren Jen Sun, 
1969, Eqs.16 and 17). 
The internal pressure P may be expressed as 

P= P, - prd os ae) 
where P, is the injection pressure and pYd is the 
overburden pressure, P is the density of overburden and Y is 
gravitational acceleration. The fracture model with the 
assumption of uniform pressurization has a stress 
concentration of infinity at the edge of fracture. When the 
fact that the stress must be finite is taken into 
consideration (cf. Eq.2.1.5), the relationship between the 
overburden pressure, the injection pressure and the internal 
pressure can be written 


2 = 
p= S0-v >) 1 (Bofrd $1} p-prd) (3.1.9) 
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where P is the average of the internal pressure over the 
entire surface of fracture. As long as the average pressure 
is assumed to act uniformly on the entire surface of 
fracture, the analytical solutions, Eq.3.1.1-Eq.3.1.7, can 
be used for the calculation of the surface deformations for 
‘this model. | 

Green (1949) and Sneddon (1946) studied the deformation 
field due to a penny-shaped fracture normal to the surface 
Subjected to the uniform internal pressure in an infinite 
elastic medium. In order to obtain the solution for the 
semi-infinite medium, we put an image source of the actual 
fracture at a symmetric position with respect to the 
free-surface, and obtain a stress field which satisfies one 
of the boundary conditions in terms of the stresses on the 
free-surface, i.e., Oy, 74=0 where n and t denote the normal 
and tangential directions of the free-surface. In this case, 
the other boundary condition, nfe Sra.Onune.O on the 
free-surface remains unsatisfied. Subsequently we have to 
solve the so called "Boussinesq problem". Namely we need to 
obtain a new stress function which satisfies the boundary 
conditions of On+~=0 and the normal load as equal and 
opposite to the normal stress resulting from the 
superposition of the original and image fields on the 
free-surface. Bo a3... 1 ét<Bqo3. hat are obtained by this 
procedure. Turning to the boundary conditions of the 
fracture plane, the image source and the new stress function 


give rise to additional stresses which are neither constant 
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nor vanish on the entire fracture plane. This implies a 
discrepancy from the initially specified magnitude and 
uniformity of the internal pressure of the fracture. 
Accordingly Eq.3.1.1-Eq.3.1.7 are not exact solutions to the 
surface deformations due to the fracture normal to the 
Surface subjected to uniform internal pressure in the 
Semi-infinite elastic medium. In order to obtain more exact 
solutions, the procedure - derive another stress function 


which reduces the additional stresses on the fracture plane, 


and solve the subsequent "Boussinesq problem" - should be 
repeated. 
Eee diet leat is an example of theoretical surface 


deformations due to a horizontal penny-shaped fracture. 
General features of the surface deformations are; 

(1) The maximum value of uplift occurs above the center of 
fracture, and the amplitude decreases monotonically to zero 
at infinite distance from the center of fracture. 

(2) The radial displacement has zero value just above the 
center of fracture and increases its amplitude to reach the 
maximum value at some distance from the center of fracture, 
and after that, decreases to zero at infinity. 

(3) The tilt field shows a pattern similar to that of the 
radial displacement except for an opposite sign. 

(4) The horizontal strain has the positive maximum value 
just above the center of fracture and has the negative 
minimum value at some distance from the center of fracture, 


and after that, increases its amplitude gradually to zero at 
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infinity. 

Figs.3.1.2a-2d show the maximum and minimum values of 
the uplift, horizontal displacement, tilt and horizontal 
Strain due to a fracture with various values of radius a and 
depth d. The uplift and horizontal displacement fields can 

be expressed in the general form 
Pearoveard ., </d } 
and the tilt and horizontal strain fields, in the form 
Pejieasa, r/d.) 
where £ depends only on the internal pressure and the 
elastic constants of material, and g denotes a functional 
relation of each deformation field and the geometrical 
parameters of a fracture and observation points. The 
function g relates to the shape of the surface deformation 
fields, and the quantity f relates to the amplitude of the 
fields. Therefore the quantity f will be designated as the 
"amplitude factor” of fracture. The amplitude factor f£ and 
the function g are non-dimensional. In Fig.3.1.2, the 
variation of a and d is given by the parameter a/d. The 
values in these figures are calculated by setting - 
feet. 0 : d= 1.0 km 
For fractures at arbitrary depth d, the vertical axis of 
Fig.3.1.2a (uplift) and Fig.3.1.2b (horizontal displacement) 
should be multiplied by the factor d (in km). 
The following descriptions hold true for fractures for 
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Pie Cel se on ees The amplitude of extrema 
of the surface deformation due to a 
horizontal penny-shaped fracture. 


Maximum of uplift 

Maximum of horizontal displacement 
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Max. of radial displacement 





= 0.2020.01 (3.1.10) 
Max. of uplift 


Max. of strain (pstrain) 


$$$ = 1.02 20.05) ~ 0.470821 .0) (301011 
Max. of tilt (pradians) * : iss 


Maxe of uplift (mm) 


ooo (1.020.2)0 (eters 
Pexemel itl te. (radians) ; 


Min. of strain 

=0.28(5=0.05) ~0.48(2=1.0) (3.1.13) 
Max. of strain 
In Eq.3.1.12, d must have the unit of km. 

These relations show that the ratios depend only weakly 
if at all on a/d. This suggests that it is very hard to 
estimate a/d (or a when d is known) by the ratios of extrema 
of different kinds of data. Eq.3.1.10 indicates that the 
measurement accuracy of the radial displacements should be 
at least four times of that of the uplifts, from the aspect 
of detectability. Similarly the measurements of horizontal 
Strain should be about twice as accurate as the tilt 
measurements, as #showNeetNurnEQ«321414. gthe exelations,ocf 
Eq.3.1.12 indicates that, when the depth of fractures is 
about 1 km, the magnitude of the tilts in microradians and 
the uplifts in millimeters is comparable, but the former 
becomes larger than the latter by the factor d (in km) for 
Shallower fractures. However the above remarks hold true 
only when the observations are near the extrema of each 
deformation field. As shown later, the locations of extrema 
of the tilt and horizontal strain fields depend on the 


fracture parameters. Therefore, in the situation where the 
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observations are restricted within a narrow region above the 
center of fracture, the uplift and the strain responses 
might yield greater detectability than others. 

Fig.3.1.3 shows where the extrema in tilt, horizontal 
displacement (or zero strain) and horizontal strain occur in 
terms of various values of a/d. We might expect that the 
location of an extremum indicates the fracture parameter 
a/d, or the radius a of fracture when the depth d is known. 
mrdes.1-3 shows that fractures with a/d<0.5 do not -.yield 
remarkably different locations, but those with a/d>0O.5 
yield some difference in locations which depends linearly on 
the parameter a/d. However, larger values of a/d do not give 
sharp extrema. Therefore, from the aspect of practical 
application, it seems hard to estimate the radius of 
fracture ( or a/d ) by finding the extrema or nodes. 

Table 3.1.1 shows estimates of the magnitude of 
Bmp cuee mactor £ for the -fractures of depth -50 -- 1000 m . 
Mie calcutation is based on Eq.3.1.9. This table contains 
three cases of AP=1 MPa, AP=10 MPa and AP=infinity, where 

AP is the difference between the injection pressure P, and 
the overburden pressure prd, i.e., 
AP = P, - ptd (3641318) 
The upper limit of f£, which corresponds to the case AP = oo 
MPa, can be estimated by the following equations. Letting x 
be the ratio P,/pYd, Eq.3.1.9 can be expressed in the form 
f= F,- fy Fs 


where 
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PAC CS ol teu ees The location of extrema of 
the surface deformation due to a 
horizontal penny-shaped fracture. 
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Table =S8h5eC¢E 644 Estimates of the 
amp |i tude factor of the surface 
deformation due to a horizontal 


oenny-shaped fracture, based on Eq.3.1.9. 
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We observe that 
a ee 
Then the upper limit of the amplitude factor is given by 


iar. - £; 


Prnax 
According to the fracture model from which Eq.3.1.9 derives, 


the radius of the fracture is expressed as a function of the 
internal pressure, as shown Eq.2.1.6. An infinite internal 
pressure gives a finite amplitude factor, but simultaneously 
makes the radius of fracture infinite. Therefore, the finite 
value of the amplitude factor does not imply finite 
amplitude of ground responses. 

In Table 3.1.1, the density Pp, Poisson's ratio VY and 
Young's modulus E are so chosen that 

mi etamenscc rani t= 10). cork =i t0'sMpa 

For most rocks, Young's modulus E ranges from 10% MPa 
(soft rock}ritol l0rteMPa.icrAs ‘tshownielbyit iko: 366. 1/Gattthe 
uncertainty of E affects estimates of the amplitude factor 
(hence the magnitude of surface deformations), whereas the 
variation of pv does not remarkably change estimates of the 


amplitude factor. 
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3.2 An Inclined Fracture 

The surface deformation due to an inclined 
(penny-shaped Or rectangular) plane fracture (Mode I, 
tensile) is discussed from the aspect of the computation as 
well as the variability of the displacement field in terms 
of the change of source parameters. 

The displacement field was calculated by numerical 
integration of the fundamental equations (Maruyama, 1964, 
page 330) for the surface displacements due to a point 
source. These equations derive from elastic theory of 
dislocation, and assume that both Lame's constants have 
equal values. The following four problems are examined; (1) 
Variations due to spacing of points by which a fracture 
plane is subdivided for numerical integration. (23 
Dependence of the surface displacement on the _ size of 
fracture plane. (3) Variation of the surface displacement 
with the dip-angle of fracture. (4) Difference of the 
surface displacements resulting from different dislocation 
distributions. , 

For numerical integration, the selection of lattice 
size is one of the problems. Table 3.2.1 shows the surface 
displacements calculated by using several lattice sizes. 
Fig.3.2.1a shows the geometry of the rectangular fracture 
( a=b=0.4 km ) used for the calculations which were made at 
places along x-axis and y-axis. Supposing that the 
computations resulting from the finest spacing ( 21x21 ) are 


nearly equal to the exact values of the displacements, the 
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crudest spacing ( 3x3 ) in this analysis yields a maximum 
discrepancy of-less than 10 percent of the exact values, and 
a spacing of ( 5x5 ) reduces the maximum aisectpetae toa 
few percent. The assumption of a point source ( 1x1 ) 
results in the maximum discrepancy of about 50 percent. The 
Spacings of ( 3x3 ) and ( 5x5 ) are about 0.26 x depth and 
0.16 x depth of the fracture, respectively. In the following 
calculations, all the grids are square and grid intervals 
are taken to be 0.08 km. 

Figwd~2.2--shows—+the-—variation—of—the-—shape ‘of? the 
displacement fields (uplift and x-componenent of horizontal 
displacement) in terms of the change of one of the size 
parameters of the fracture. The calculations are based on 
rectangular fractures (horizontal and vertical) with a fixed 
value of a (=0.4 km) and two different values of b (=0.24 km 
and 1.68 km). The profiles in Fig.3.2.2 are along x-axis in 
Bigeaec. 

Fig.3.2.2a shows the normalized amplitudes relative to 
the maximum value for the uplifts due to a horizontal 
fracture. The horizontal displacements due to a horizontal 
Or a vertical fracture and the uplifts due to ae vertical 
fracture respectively do not show maximum amplitudes at 
fixed coordinates, which do not necessarily coincide with 
one of the calculation points specified. As a reference, the 
horizontal displacements due to the horizontal fracture, 
normalized by the maximum amplitude of the calculations, are 


shown in Fig.3.2.2b. Fig.3.2.2c shows the uplifts due to the 
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vertical fracture, normalized by the amplitude at the place 
above the center of the fracture. 

It can been seen from Fig.3.2.2 that, without regard to 
the difference of the size parameter b, the resulting 
displacement fields have an almost similar shape along the 
profile, which is perpendicular to the direction of the size 
Parameter b. In order to examine the similarity of the 
shapes of the two displacement fields resulting from 
different values of b, we may evaluate the cross correlation 
of these two data arrays. The cross correlation for each Be 
the four different displacement fields is shown in Table 
3.2.2, which demonstrates that the shape of the surface 
displacements along the profile depends only weakly if at 
all on the size b. It is also observed from Table 3.2.2 that 
the vertical fracture can reflect the change of the size 
Parameter b in the shape of the surface displacement fields 
more sensitively than the horizontal fracture Hoeeh. 

The above observations suggest that, as long as 
observations are made along such a profile, they do not 
provide enough information to determine the size parameter 
b. In order to estimate this dimension of the fracture, 
observations in the direction of fracture strike (y-axis in 
Fig.3.2.1) should be available. On the other hand, the above 
mentioned weak dependence implies that a 2-D fracture model 
can be used for the determination of. the other source 
parameters. Even if the overall shape of a deformation field 


depends to some extent on the variation of b (for example, 
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Displacement 





ety Vines ates Re. co, Cross-correlations of 
two displacement fields resulting from 
different values of the size parameter b. 
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the uplifts due to a vertical fracture shown in Fig.3.2.2c, 
which yield the lowest correlation in Table 3.2.2 ), the 
Shape of the deformation field over a relatively small 
measurement area (for example, x=0.0 - 0.2 km) can become 
independent of variation of b, and the relative amplitudes 
can be expected to be conserved. If the deformation field 
over the small area responds well to changes of the other 
source parameters such as length a and dip-angle of the 
fracture, 2-D fracture models can be used for the 
determination of these parameters. It should, however, be 
taken into consideration that the small size of a must not 
result in significant change in the shape of displacement 
fields due to changes of its value, and that the estimate of 
displacement dislocation or driving pressure directly 
depends on the absolute magnitude of observed displacements. 
The maximum amplitude of the displacement fields due to each 
of the fractures with various values of b is shown in Table 
ad. Jt 

Fig.3.2.3 shows the uplift, horizontal displacement, 
tilt and horizontal strain fields due to a penny-shaped 
fracture with a dip-angle of 60 degrees, radius 0.2 km and 
the geometry shown in Fig.3.2.1b. The latter two fields were 
calculated by applying bicubic spline analysis to the former 
two fields, respectively (IMS Library, Subs. IBCICU and 
DBCEVU). : 

Figs.3.2.4a and 4b illustrate the profiles of the 


uplift and the horizontal displacement (x-component) along 
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Fig.3.2.3 ..... Surface deformation fields 
due to an inclined penny-shaped fracture 
with a dip-angle 60 degrees and 
characterized by a constant displacement 
dislocation. 
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UPLIFT 
DUE TO AN INCLINED PENNY-SHAPED FRACTURE 





-1.0 -0.6 0.0 0.56 1.0 KM 
X 
Bi «ta eatin be» « fae Variation of surface 


uplifts due to an inclined penny-shaped 
fracture with various dip-angles. 





HORIZONTAL DISPLACEMENT (X-DIRECTION) 
DUE TO AN INCLINED PENNY-SHAPED FRACTURE 





-1.0 -0.5 0.0 0.6 1.0 KM 
X 

Fig.3.2.4b hese. Variation of surface 

displacements (horizontal) due to — an 

inclined penny - shaped fracture with 


various dip-angles. 
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the x-axis due to a fracture with the same geometry as the 
one shown in Fig.3.2.1b and with various dip-angles (Dip=0, 
30, 60 and 90 degrees). These figures show that the 
displacement fields due to an inclined fracture are very 
sensitive to the change of dip-angle. 

The assumption of constant displacement dislocation on 
the plane of fracture is at variance with real conditions. 
Here the displacement dislocation will be given by the 
Symmetric and linear distribution 

Deere Dea ADiri +0-s<--lrl spa 
where r is distance from the center of fracture to the 
center of grid cell, and a is the radius of fracture 
Pere 8.2 te) 

Fig.3.2.5 shows the difference between the magnitudes 
of the surface displacements due to each of the two model 
fractures, one of which is characterized by the constant 
displacement dislocation and the other is by the symmetric, 
linear displacement dislocation. The geometry of the former 
model fracture iS same as one in Fig.3.2.1b, and the 
distribution of displacement dislocation for the latter 
model was given by D, =30 mm and AD=30mm/0.2km, which 
yields an identical cavity for the svoemoded fractures. 

In Fig.3.2.5, a- Square in each grid represents a 
calculation point, and the shadow shows the degree of the 
difference in percent between the amplitudes of the 
resulting two displacement fields. A large difference of the 


amplitudes occurs at the places where the displacement is 
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very small or the spatial gradient of the fields is steep. 

In general, changes in the displacement dislocation on 
the fracture plane do not significantly change the overall 
Shape or the amplitude of the surface displacements. This 
implies that appropriate selection of the observation points 
for an inversion analysis permits us to adopt the simple 
fracture model characterized by constant displacement 
dislocation. It is also seen from Fig.3.2.5 that, if strike 
and dip direction of a fracture can be inferred roughly in 
an experimental site, observations over the half-area in the 
dip direction across the strike above the fracture center 
Meenmetne region of x= -0 in Fig.3.2.1) might provide 
useful information for an analysis based on the simple 
dislocation model. The assumption of uniform pressurization 
in a fracture results in the separation of fracture walls 
having elliptical cross section. The characterization of 
displacement dislocation distribution by an elliptical form 
can be easily assumed to yield displacement fields similar 
to those resulting from a constant displacement dislocation 
distribution. 

The above results are concerned with a particular case, 
(a typical length of a fracture)/(depth)=0.4. The 
variability of the surface deformation fields in terms of 
changes of the source parameters such as dip-angle and 
dislocation distribution should be examined in view of the 


methods which will be used for inversion analyses. 
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3.3 A Shear Fault Model 

In the previous two sections, surface deformation has 
been discussed by assuming that a tensile fracture will be 
produced by hydraulic fracturing. As described in Chapter 2, 
there is also a possibility that hydraulic fracture breaks 
out as shear failure or as combination of shear and tensile 
failures. It is necessary to understand the ground response 
associated with the presence of shear cracks. 

Analytical expressions for the surface displacements 
due to an inclined rectangular shear fault model with 
constant slip have been obtained by Manshinha and Smylie 
(1971), under the simplification that both Lame's constants 
are equal. An example of the surface deformation fields due 
to this fracture model is shown in Fig.3.3.1. The tilt and 
horizontal Strain fields were calculated by applying bicubic 
Spline analysis to the uplift and horizontal displacement 
fields, respectively. 

We have to observe that the tilt and uplift fields in 
Fig.3.3.1 show roughly circular shapes. If we assume that 
hydraulic fracturing must always cause tensile failure, we 
would incorrectly infer from such roughly circular 
deformation fields that the induced fracture must be 
horizontal and have a circular or square shape. The example 
in Fig.3.3.1 also demonstrates that measurements of the 
horizontal displacement and strain fields are indispensable 
to determine source parameters of the fracture. These 


Suggest that use of a particular kind of deformation 
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FRULT PARAMETERS 
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Fig.3.3.1 ..... Surface deformation fields 
due to an inclined rectangular shear fault 
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measurements, limitation of the number and inadequate 
deployment of measurement points can not provide enough 


information to delineate the form of fractures. 
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4. Determination of Source Parameters of a Horizontal 


Penny-shaped Fracture 


4.1 Basic Procedure 

-One of the possible measurements at a measuring point 
at distance r from the center of a fracture of radius a, 
amplitude factor £ and at depth d can be expressed in the 
form 

Garam ae gd +) ees) 

As seen from Eq.3.1.1-Eq.3.1.7, the function g (and G) are 
highly non-linear with respect to the unknown parameters, a 
and d. The depth of fracture, d, may be usually regarded as 
a known parameter because the depth of formation of the 
fracture is controlled. However, there is no evidence that 
the fracture is induced at the position centering around the 
depth of fracturing operation. Therefore it is more general 
to regard the depth of fracture as one of the unknown 
parameters. 

For simplicity, it is assumed that the depth is a known 
Parameter. The function g is still non-linear with respect 
to the parameter a. Here we observe that the parameter f 
appears linearly in the observation equation. We assume 
several arbitrary values of a, for example, at an increment 
of Aa, and apply the least-squares procedure to a set of 
measurements. We can determine f for each of the assumed 


values of a. Then we obtain the least-squares errors 
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resulting as a function of a, E(a). We may expect a 
particular value of a which yields the minimum value of 
E(a). This value of a might possibly indicate the rough 
estimate of the actual radius of the fracture. If the 
observed data is noise-free, this particular value of a 
should indicate the exact value ( within the certainity of A 
a ) of the fracture radius. In practice, we can not expect 
any field data without noise. Here I discuss the degree of 
resolution for the fracture radius when measurements are 
disturbed by noise. In other words, the problem is to 
understand how well the surface deformation delineates the 


changes of the source parameters. 


4,2 Resolution - Basic Concept 

Let a, and £ respectively be the true value of the 
radius and the amplitude factor of a fracture considered. 
Suppose one of the possible measurements at the j'th 
measuring point of distance r yields the value G' due to 


J 


the fracture. We can express G' in terms of the actual 


° 9 
response of the surface, 93. , and a "noise component", N, 


J 
in the form 

The = ets N. Ea 
where 
Gawd ae; rj ) 


When we assume an arbitrary value of a, we have the model 


equation 
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in which the unknown parameter f appears linearly. Then we 
can apply the least-squares procedure to a set of 


measurements G'; , j=1~J, in terms of the model equation. 


E¢asja = SCG Hg )*_S°6 
a 


i ea ie 
oE/da= 0 
It is straightforward to show 
— fae y2 rey Y (g 2 e e 
j 3; 


Substitution of Eq.4.2.4 into BG. 4.2.3 yields the 
corresponding least-Squares error 
FG? 946) e 2 
E(a,a) = >, { G' - —— 5, } (4.2.5) 
j ¥. (4) 
5h a0 e i 
£5, +N) 3. | 


(136) 
e 2 
De C3) 


00 2 
= 2, (F9,+Ny) 


Using the notation for inner product 


<A, B> = 2,048; ) (4.2.7) 
Eq.4.2.6 becomes 
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Substitution of g = g into Eq.4.2.8. yields 
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<3,N> 


Etats y= CNN 
<3,3> 


(4.2.9) 


E(a.,a.) corresponds to the least-squares error with respect 
to the true value of the fracture radius. From Eq.4.2.8 and 


Eq.4.2.9, we have 


fs ve °o Lh 
Bta,.a) - Ela, 53. )}4= ag <9. 9)— eras ) 
{853 
és . 3.95<3,N > (4.2.10) 
+2 0 Lj ot EES 
as te 


In Eq.4.2.10, the first term on the right hand side is 
always positive, whereas the sign of the second term can not 
be determined definitely because this term depends on the 
Mature of function g and noise N. This suggests that E(a,,a) 
does not necessarily have the minimum value at a=a, (a, is 
the true value of the fracture radius). 

Fig.4.2.1 schematically illustrates E(a,,a) and the 
concept of the resolution relating to this inversion scheme. 
In the practical analysis of the actual measurements G', the 
curve of E(a.,a) is obtained by calculating Eq.4.2.5, and 
eventually represented as discrete data series corresponding 
to discretely sceceee values of a. 

Suppose we know the noise level, <N,N> . We may 
represent the noise level by a horizontal straight line as 
shown in Fig.4.2.1. The curve of E(a,,a) and the line <N,N> 
should yield intersections, for example, at a=a, and a=a,. 


We can not judge which point on the curve of E(a.,a) 
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corresponds to the true value of a, but from Eq.4.2.9 we 
understand just that the true value of a (=a,) must lie so 
that 
So ae ee eS 
If the noise level is high, for example, as illustrated by 
the broken line in Fig.4.2.1, 
Ores 7a 9"S fa; 
In this case, the lower limit of the estimate of the 
fracture radius equals zero. 
Here we will examine the range of a which satisfies the 
condition 
Bay fae ao as.” 4027113) 
Eq.4.2.6 indicates that E(a.,a) has the minimum value zero 


when 
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where & is introduced as a constant. Substituting Eq.4.2.12 


in the form 


g, = -N,/ 4 or § = CEAG oN ao (4.2.13) 
into the second brace of Eq.4.2.8 , it follows that 
£3,385 
° oO 5 
EC diya auc bantcg. 9% ome 
S547 (4225745 
<3,N> 
~ (<N,N>~ SSHEZ ) =E’ 
€ 353 


E' represents the minimum of E( a.,a ), which results from 
Eq.4.2.12. In general, the coupling between the noise and 


Signals shown in Eq.4.2.12 may be excluded from the 
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consideration. Accordingly, the inequality 
EC a pat.) 2,465" ier e2e15) 
is always satisfied. From Eq.4.2.9 and Eq.4.2.14, we can 


express the condition to satisfy Eq.4.2.11 in the form 





atar> e Zz 

ose % {323 <IoN 

£7(<8,§)- ae )- sed her yg < 0 (ae) 
539 £353> 


Neglect the term of <g,N> 2/<£9,9> (this is positive value), 
and rewrite Eq.4.2.16 as 


£3,9> 
Oo oo 3 
f*( <3,3>- ee 


<39,3? 
It is obvious from the above procedure that the values of a 


) < 2<N,N> (4.2.17) 


which satisfy Eq.4.2.17 simultaneously establish Eq.4.2.16 


and Eq.4.2.11 as well. Let 


etre 2 ’ 
{ o °9 Sy 2 
aa Ae ey. ale (4.2.18) 
3 hint bee 4 KGaeh.S 
Oo - ~< N,N>| : | (4.2.19) 


Eq.4.2.17 can be expressed in the form 


6 
fF, Gas,ay) eto, (4.2.20) 


Oqg(Ge,d) can be regarded as an indicator how the 
function g=g(a,r) changes its overall shape (g; ,j=1-J) in 
terms of the changes of a. Og (aed) represents the degree 
of similarity between the figures of g=g(ae,r; ) and g=9(a Tj 
), j=1-d. For example, let J=2, and 


. ° ° e 
g, = g(a ,r, ) ,= glataa,r, ) 49, =.9,..5.9, 
e ° ° e 
g = g(a ,r, ) g,= glataa,r, ) Aes. 
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Then Og is given as 


Ag, - 9,49, 
26; = S45. 3.43 
g° + 64, 


Consider the two figures shown in Fig.4.2.2. It is clear 
that the two polygons, ABCD and EFGH, are similar figures. 
Indeed, substitution of the values indicated in Fig.4.2.2 
yields S =0. If J=1, the "Shape" of g is just a line, and 
regardless of the variation of a, the shapes of g are always 


similar figures. Hence 9, =0. Eq.4.2.18 will be expressed in 


the form 
oO @ 9 
Be <9,9> 
TERS @ahee: / UNS eres — | 
£9,9>< 9,9) 


The second term in the brace is related to the 
cross-correlation between the two data arrays, g and g. 
Accordingly Gg can be also considered as a sort of 
correlation function of g in terms of the dual variables, a 


ana xr. 


g(a,r) g(a+aa,r) 
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4.3 The Resolution Related to the Model Functions g 

I will evaluate Gj using the model function g for the 
surface deformation fields due to the horizontal 
penny-shaped fracture, and examine whether or not the 
function g has the sufficient characteristics to satisfy the 
condition expressed by Eq.4.2.20 under a given noise level. 
It 1s noted that the theoretical estimate of the resolution 
should be made with respect to /26,, whereas the resolution 
( error estimate ) for the analysis of actual observed data 
should be evaluated against OQ). 

Eq.4.2.20 indicates that large a variability of 0 
with changes of a is preferable in order to determine the 
radius of fracture accurately. Og(ao,a) depends on the 
fracture parameters, a and d, and the geometry of measuring 
points, ri j=1~ J. We need to understand how the function 
93 depends on these parameters. We are free to design the 
Site configuration of measuring points, if the situation 
permits. It iS natural to choose such a configuration of 
measuring points that yields as high a resolution as 
possible in the determination of fracture parameters. 

The analytical evaluation of 3 is difficult. In order 
to obtain Og(a.,a,J), the model field of g(a,rj), jni~ J, 
was generated by the following specification. 

f = 1.0 d = 1.0 (km) 
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The interval of measuring points (Ar) = d/50 
and ye Oe es oe 
In the following, the variables a, and a of the function Og 


(a,,a,J) appear in the form 


Og( a, ,a,+faa,J ) 
ao=ald , Aa=d/20 , R=0, 41, +2,-- 

The function 4, for. each of the uplift, the tilt, the 
horizontal displacement and the strain is identified by 
aeeo, , Op and 6, . 

Since the calculation of Og is based on the model data 
generated by f=1.0 and d=1.0 km, comparisons of Og for 
arbitrary values of f£ and d with noise level,(Eq.4.2.20), 


are made by 


Fade Sin FLOR, 
rd uO, ween 12) Oy 
Matas, “aie Ci 302) 
£07 < 2 Oy, 
Oe mete On 2 
where G,,, (g=U,T/R,S) denotes the noise level for the 


respective deformation field. 


Rag.4 <3 yl shows O7(a.=0.4,a,J=50). This figure 
illustrates general features of most Og with fixed values of 
a and J. The following features can be seen. 

(1) Gg behaves almost linearly with respect to tha, 
Sees ae in the region of 2> 0 and small values of -Q 
Gg «x + haa 


(2) In general, Og has steeper slope for 0>0 than that for 
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OG 64s, ay-£48,J) < O5(a,,a,+24a, 7) 


" " 


This suggests higher resolution of the upper limit of "a 
than that of the lower limit. 

Fig.4.3.2 illustrates the dependence of Og on J (number 
of measurement points). Og has the steepest slope for a 


Particular number of J, ( <Ji Jee The following relations are 


found. 
50) oe 7.0 for G, 
a 50 ~ 80 For OF 
Paelee WC Rete Wc 
60 ~ 90 for On 
pee 5004 SF oreo. 


where a/d = 0.2 ~ 1.0. Since the parameter J corresponds to 
Span length of measuring points, the above relations can be 


replaced by the relations between the length of the span, L, 


and the depth of fracture, d. 


(A804 Fa for O;, 


hd x: a ( 1.0-1.6 )d for O& 
Peas | XAT = Ik (4.3.4) 

30 Ceteee 1 OO. Mears 

202.0~1 097d BErur Oe 


As shown by Eq.4.2.20, the larger value of 6g is more 
desirable in order to obtain high resolution. Therefore x 
will be referred to as the "optimum measuring span". It is 
obvious from the nature of this inversion method that a 
longer span of measuring points does not always yield higher 


resolution in the determination of the fracture parameters. 


The magnitude of the surface deformations all decrease 
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toward zero at large distances from the center of fracture, 
and the overall shape of the surface deformation field 
becomes more "Similar". 

The above relations, Eq.4.3.3 or Eq.4.3.4, are 
illustrated by Figs.4.3.3a-d. These figures also indicate 
that O5(ao,ao-Aa,J) of larger a, changes its value more 
drastically than that of smaller a, with changes of J. This 
feature is significant for O;and OR. 

tie ratio of Oglae,ao- Aad) / Og(a.,ao-Aa,J=10) is 
tabulated in Table 4.3.1. Suppose the tilt responses 
associated with a fracture at depth 500m and of radius a 
=100m ( a,=0.2 ), and the noise level such that the 
resolution of mp ete, at F190; 4.5.48) for the ‘fracture 
radius is attained by the optimum span of measuring points ( 
L=500m ). Then Table 4.3.1 ( see * ) indicates that, if the 
span length is J=10 ( L=100m ), the attainable resolution ( 


Be in Fig.4.3.4 ) is as low-as 20.9 m. 


Next, comparisons between G7, ay iy ee O. wre 
shown. The relation of Eq.3.1.12 has shown that the signal 
levels of uplift and tilt fields in terms of maximum 
amplitude are related by the factor (1.0+0.2)d in the ratio 
of uplift/tilt. On the contrary, the analysis of resolution 
factors results in the relations 

d Oyla, »a,-Aa,J=10) 


= 2.00 (a,=0.2)~4.98 (a,=1.0) 
O7(a, »8,-Aa,J=10) (4.3.5a) 
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dO(a,,a,-da,d=d ) 
a = 0.36 (8.=0.2) ~ 0.52 (8_=1.0) 

O7(a,,a,-Aa,J=J ) % ‘ (4,3.5b) 
where d has units of km, and the ratios have units of 
mm/microradian. 

Assume that the noise level of the uplift data in 
millimeters is comparable to that of the tilt data in 
microradians. Eq.4.3.5b suggests that, if measuring points 
are arranged to the optimum span, the analysis of the tilt 
responses can yield a more unambiguous solution than that of 
payee upitet data, as ~‘long “as * the™ depth” of fracture’ 1's 
Shallower than 2 km - 3 km and the ratio of radius/depth is 
less than 1.0. On the other hand, Eq.4.3.5a suggests that, 
if measuring points are limited to a narrow area above the 
center of fracture, the uplift data can be more useful than 
the tilt data for the analysis of fractures at depths deeper 
Pane ooo m — SOD i. 

The ratios of O; /O7 show the relations 


O3(a, a ,-da,J=10) 


ed dae (a,=0.2)~ 7.94 (aj=1.0) 
G7(a,,a,-Aa,J=10) 

ra (4.3.6) 
O5(a, .a,7Aa Jad ) 


) 


O26 (a,=0-2) ~0.71 (ay=! 0) 


Fai} 


O7(a, ,a,-4a,J= 


Oo 


where the units are microradian/microstrain. The ratios of 
OG, /O, show the relations 


O(a, a -4a,J=10) 
eee Ve = L556 (a,=0.2) ~28.2 (a,=1.0) 
Gala, 1457 4a, J=10) 

A ale Sy 


Gi(a_,a.-aa,J=d ) 
ie ee eae 0, 2) ~ 2. 12 a0) 
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where the units are mm/mm. The above relations all show 
different feature from the relations of Eq.3.1.10-Eq.3.1.12. 

Several examples of quantitative details of Og are 
Shown in Figs.4.3.5a-5h. Table 4.3.2 particularly tabulates 
: Og (a.=0.4,a.-Aa,J). Assume we have surface deformation 
fields due to a fracture of radius 160m, amplitude factor 
10-° and at depth 400m. Table 4.3.3 gives the upper limit of 
the noise level which will allow determination of the radius 
of fracture within the certainty of +20m from the 
Measurements ( Ar = d/50 = 8m ). The values in this table 


are calculated by using Table 4.3.2 and Eq.4.3.2. 


So far, the resolution functions Oy have been evaluated 
by setting a configuration of measurement points, which 
Should be discussed as a variable of O,. Here supposing the 
five types of measurement configuration shown in Fig.4.3.6, 
I examine which configuration yields the maximum resolution 
for each of the deformation fields of uplift, | horizontal 
displacement and tilt. The results are summarized in Table 
4.3.4. The optimum configuration, which yields the highest 
resolution, 1s shown in the column of "Best". Assuming that 
a resolution 20 percent lower than the highest resolution is 
the minimum acceptable measure, the corresponding reasonable 
configuration is shown in the column of "20 percent 
allowance". 

The comparison of the resolution factors attained by 


the respective optimum configuration for the three 
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limit of noise-level «to attain a 
resolution of +20 meters for the radius of 
fracture. 
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span length of measurement points ( depth7! ) 
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acceptable measurement configurations. 








a 
= —— ar~rsipneaieiteaeiiagaeaciemasiaae 
eeatinemasysiiaiantoats ———————— 
; ; : 
‘ eihivieae orn 
Coté > os T2368 x 
; 
r — 
i - - Oot . 
: 
‘ ‘ - ! ' ¢ wa it & = 
i i > J 
§ Py t 
t + | ¢ = t 
| 
ealiincea a canals rr 
| 
: 
: ivy ~ j~ 3@ i 
: : 
: él | = | } 
; 12 
——— - Ot cg nn nee lh ne 
iy ; ' 
} 
: n> o~<6e M 
¢ a < 
; i ¢ 
-$9¥ 
~~ AD \ 
* ¢ | 
’ a 1 
ps | oa. 
e [eracsiieoh 
" ) a a a: — : a 
| inomeasciqeid . 
a es : 
ot > bX | 1 2 ose Seo. | a 
me ' ' i 
| | 


| ‘isi oft ze omee bEi~Qiit = J 










oa f 
= 
ry a 
a oe Lh. ont ee Ee te te ee 
oy 


) #:0ioR1t- Yo ewiber 
7 ' - AV, hae 8 “piwtaerd to diqsh Bp ob 
Ra =) voces “Yo 4 Nt 
”.] : 





140 


deformation fields indicates 


Gr Gu 


Or Or 
her tne fractures which 0.2 < a/d < 1.0. This analysis 
also indicates that resolution factors for the tilt and 
horizontal displacement fields depend heavily on the type of 
measurement configuration, especially for the fractures of 
large a/d. This probably results from the fact that the 
location of the extrema for these fields shifts remarkably 


with changes of the radius. 


In the analysis of observed surface deformation data, 
the amplitude factor £ is obtained as a function of the 
variable a (see Eq.4.2.4), like the least-squares error E(a) 
(see Eq.4.2.5). If we can estimate the order of magnitude of 
the amplitude factor of actual fracture, we may expect that 
the behavior of f(a) can determine a possible range of the 
fracture radius. The behavior of f(a) may be evaluated by 
the calculation of 
3, 9) 
3, §> 
Similar to Og(ao,a,J), f will be expressed as Fg 
(Ae;ayd). 

Fig.4.3.7 shows ae aida) ae The values in this 
figure are calculated by using the model tilt field 
generated by the previous specification (Eq.4.3.1). The 


6 
following features of fg /f are found. 
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Pig. 4 et eee ee The amplitude factor as a 
function of an unknown parameter a, 
normalized the true value of the amplitude 
facture, 
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Mirueres / f of smaller a, has steeper slope. 

(2) In the region of 2<0 ( angers: ares of a ), it shows 
more distinct difference from the actual value of the 
amplitude factor. 

(3) The behavior of eye does not depend on the kind of 
data (g=U,T,R,S) and the number of data points (J=10~ 
j3=100) and all show an almost same curve as Fig.4.3.7. 

(1) and (2) suggest that, if the order of magnitude of 
tne actual amplitude factor can be estimated, the resulting 
values of f(a) determine the lower limit of possible values 
of the fracture radius, especially for fractures of small 
radius compared with the depth of fracture. It is noted that 
the dependency of fa PE on a, and +2aa shows the opposite 


feature to that of Og - 
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4.4 S/N Enhancement by Combination of Two Kinds of Data 

There are several measurable quantities, such as 
uplift, tilt, radial displacement and strain, relating to a 
fracture. When some combination of these measurable 
Quantities is considered as input data for the inversion 
problem, we may expect some enhancement of the signal to 
noise ratio. Here I will consider the simplest case in which 
two kinds of observed data, P' and Q', are combined in the 
form 

Z' = Pl; oy xO" 5 (4.4.1) 

where j denotes the j'th measuring point, and x is 
introduced as a variable parameter (it will be called the 
"stacking factor") and assumed to be independent with 
respect to j. Incorporating the combination of Eq.4.4.1 and 
the inversion scheme described in the previous section, I 
will obtain an equation for the determination of a specific 
value of x which will yield higher resolution than that 
achieved by the observed data P' or Q' alone. 

Consider a fracture of radius a, and amplitude factor 
E and express the observed values of data P and Q at the 


j'th measuring point of distance r;, in the form 


Ny, 
° ( ) ? ° aC 
Pl: om Serine) oe Daa ING 
J “eS Sica (4.4.2) 
° °° #698 (qy e ° 
! com pe j\ 
Q Aidt d qa (a> oP) Yes q; x3 N5 


where p and q are the model fuctions for the deformation 


fields due to the fracture, respectively relating to data P 
(P) (q) 
and Q, and N and N are noise components of data P and Q, 
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respectively. Define a function 


Z “Ailes xd 
and let 
9 °o ° 
i has i alam (404.3) 


In the same way as in the previous section, the following 


notations are used. 


B = pla r,) (Lous) 

2 = q(a,r.) (4.4.5) 

25 = 2(a,r) (he a.G) 

eee Pt AB) (4.4.7) 
ey res 


Denote GO defined by Eq.4.2.18 for the functions p,q and z 


by 6,, Og and O63. Also denote G for the noise components 


by Bet Ona and O,,, respectively. These ©@ are given by 
} ole e.h.° 
JO.* = <B,p>-<B,BY/KB, BY (4.408) 
2 oe) 
JOS = <954)- 48 ,49°/4 54> (4.4.9) 
JOs= = (2,2)=42,2)°/K2 ,2> Cd ay 
(P) cP» 
JO iy = BGN aN} Ch gh31 1) 
(q) CQ) 
IOF g = <N,N> (4.4.12) 
Jor, = (hy oeax Roa x2C, 1 eee 
ie yl leetts 15) 


In the following, it is assumed that 


CP) ca) 


< N,N? 35 O ec teg 4) 
Substitution of Eq.4.4.3 and Eq.4.4.6 into Eq.4.4.10 yields 
the expression of O,.in terms of p and q, 
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where 
5. 40 e 86°88 a ze 2 
Ag = €P»P><PsP>—-<PyP> 
° ehUe ejU6©°8 Oo °9 
A, = 2( <B,4<B,P>+(B,4)<P, P> 
°o @ o oe eo 
~<P,P)<P,G>=(P,P><p,g> ) 
ie} 0° e ée fe} fo} cs . (o no) e@® @e 
As = <DsP>(4s4>+2<P54)< Dy d 7+ <4 54K Py PY 
oe e o @ e ole ° le 
= (By) =<B,d)~=24D,a)X D4) 243, DX DO 
°o O° ee e oF 8 e he 
= <B,d <9 4 = <Psd)K 450) Out. 16) 
Ay = 4,4><4,4>-<8,89° 
Bo = <D, p> 
B, = 2¢p,a) 
B, = (4,4) 
Let 
Phen oe . 
( es oy /, 207 6p CA a) 
pa OF % e es 
(RB), af f2Ori y (ive telloy 
(R= Orie , (464019) 
Beestatution of £Eq.4.4.%39/ andg)Eq.4.4.15 “into Eq.4.4.19 
yields 
1 (Asta, X+ASx“+A, x7 +A, x4 2 
(R), oa 2 = = 2 (ipetieeo) 
a oo ts Bie ee ) 
where 


_ Ce i 2 
Dy = oie 3 Ds = 4°) (4e4e21) 
(R)z expresses the resolution resulting from the combination 
z=p+xq. (R)z maybe behaves such that it has the maximum 
value at a specific value of x, which is the most desired 


"Stacking factor" for the enhancement of resolution because 
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larger values of (R)z correspond to higher resolution. 
Eq.4.4.20 indicates that, if the information about the ratio 
On p/ Sng is available, the desired value of x can be 
estimated by examining only the character of the functions p 
and q. 
The degree of the resolution enhancement resulting from 
the combination can be evaluated by the ratio 
(R) /(R), or (R)/(R), 
The ratio is understood as follows. Fig.4.4.1 illustrates 
the resolutions relating to the function p alone, q alone 
and the combined function z=p+xq. Here we assume _ the 
linearity of O-curves, which has been described in Section 
4.3. The resolution in the determination of a (radius of 
fracture) can be expressed by a range between two 
intersections of O-curve and OO, -level. The resolution 


factor b shown in Fig.4.4.1 is given by 
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Denote b for p, q and z by bp ebopand seb ea The raulG 
bz / bp is written in the form 
Pz Sh One Ont p Pp (hele 23) 
b., Der, mat wate 


and expresses the degree of the enhancement of resolution. 
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4.5 Determination of Time Varying Fracture Parameters 
Denote the radius and amplitude factor of a fracture at 
time Mexby ajiandet ft, § 
aye Sftt, 4%, wplade op§ (4.5.1) 
and let G, and 3.44 be one of the observables at the j'th 
measurement point of distance r and the corresponding model 
function, respectively, 

Deen Ghee le ny = 7,9 68 FD 1455.2) 
where, for simplicity, the depth of the fracture has. been 
assumed to be a known parameter, and noise components are 
not considered. The least-squares method with sweep 
procedure for the non-linear parameter a ( This procedure 
has been described in Section 4.2, and will be denoted by 


LSM* ) gives the redation;vfordthe datassettatitime ®t}; 


£ ms rt ee, Cee PORES 
€ 9+ Ie? 
Similar to Eq.4.5.1, let 
a, = act) ; S2.j = GC a2,4; ) 


(4.5.4) 
Peete: 8G 7 ea Pea Caan 


The LSM” gives the relation, for the data set at time t,, 


ee S< S21 G2? (4.5.5) 


< $2, 92? 
Then we can calculate the change of fracture parameters 
during the time= t,~t, as 
af* = f, - Ff, (4.5.6) 


i = a Oy (4.5.7) 


If we have several sets of observed data at various times 
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during the growth of the fracture, the fracture source 
parameters can be interpreted as time varying functions. 

Next, one of the alternative methods for the 
calculation of 48 and Af will be discussed. We have the 
“equations 


Gey = Fe ej OF) (4.5.8) 


and 


AG = Great, | = Gr, 
9G] + sade pans 5) 
An OF It7 | 


a) 
year SE stat Erle (4.5.10) 


Pa-*the ‘appropriate **initial’ value “of+ fFocsand®*'a”™ can be 
specified, and the observed data at more than two measuring 
points are available, 42 and aft can be calculated by using 
Bast 25 510. 
Define © and E by 
e, = at g, + aaf,- cde! mt AGS 


t, 
oa It. (4.5.11) 


PE NN me ee 


Apply the least-squares procedure to Eq.4.5.11 








= 0 (255.12) 
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69,99<aG, 2h >-<4G,9>¢ 9,92 > | 
Aa= (4,9><89 22) <g, 28 9 i (4.5.14) 
where numerous subscripts have been omitted. Although this 
procedure can be an alternative to the method described in 
Eq.4.5.3-Eq.4.5.7, it gives rise to several problems, which 
will be discussed in more detail. 
A first remark is concerned with the initialization of 
a and f. Eq.4.5.13 and Eq.4.5.14 require the initial values 
such that a@*0O and f+0QO . This is so because the 
calculation of "a" requires non-zero value of £ as shown in 
Eq.4.5.14. It suggests that an additional tool (such as LSM* 
) should be prepared in order to determine the initial value 
Seed eaands~t o,lt oth@s-initiale.values, are, not «specified 
correctly, the error which might accumulate should become 
serious. Usually the initial value can not be specified by a 
unique correct value but by some bounds. Therefore it is 
necessary to examine the results obtained by different 
possible initializations. 
A second remark is concerned with computational effort. 
Suppose that we first calculate 4@, and obtain the resuit 


of A@=0. In this case, we can use the equation 


af = Sebi (4.5.15) 
war Lee os 


instead of Eq.4.5.13. Eq.4.5.15 is obtained by substituting 


A@=o into Eq.4.5.10. Probably we would first compute only 
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the numerator of Eq.4.5.14 for the calculation of Aa. When 
the result is 4@ =0, we do not need- to compute the term 
< og 2955 in the denominator, and we can simply determine 
Af by Eq.4.5.15. This suggests significant reduction of 
computational effort. Similarly if we first calculate Af 


and obtain the result af =0, asaacan be calculated by the 


equation 
hee vans! | 
Aa = ae CRs (4.5.16) 
<g 2h) 


A third comment is concerned with the inherent error 
contained in the solutions, Aa and Af, obtained by 
Eq.4.5.13 and Eq.4.5.14. The true value of Af is given by 


£G2,952? © CCAS 


i ae oe ee 
at fa 425.,922 ie GT fae SAE ae 2 


(4.5.17) 


Suppose that each element in £Eq.4.5.13 and Eq.4.5.14 is 
given by quantities in terms of time { and 5. 
Substitution 

AGeel Fyizig Sq. $G7. 7,5, | Gop reret g, (4.5.16) 
into Eq.4.5.13 yields 


Apo bene tae na fos 9} 
where 


2£5,,9,>% 91,5 > 


A RR eRe RTE ree (4.5.20) 
A Sases 5,44) ta.) 


AAT = 


AAF represents the error originating from 8&q.4.5.13. 
Similarly we can imagine that 42 of Eq.4.5.14 also yields 
some difference from Aa*, However, af in the case of 


Eq.4.5.15 yields the correct solution, because the basic 
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equation, Eq.4.5.10, is a simple linear equation with 
respect to the parameter Af only. 

Here, let us Lg Pe the.-error.«source of AAdAf and 
construct a procedure to reduce the error. The error comes 
from the treatment of Eq.4.5.10. The original equation of 
Eq.4.5.10 is 

AGI = Groes = G,. 
af g,+ £49 + Af 4g (4.5.21) 


Eq.4.5.10 is an approximation of Eq.4.5.21 as the result of 


u 


regarding 
[eat ASi leis, ele ATS tet (4.5.22) 


and discarding the higher terms of 





3g sa ag : 
AG = ey eee yee (455.23) 
hence 
ag 


In order to deal with Eq.4.5.21 as a linear equation, we 
adopt the approximations Eq.4.5.22 and Eq.4.5.24. Here, it 
is assumed that only the approximation Eq.4.5.24 is 


Satisfied. Eq.4.5.21 can be rewritten into 


AG = Af ( G, +49) + 49 Ff, 
= Af-9, + aqf, (4,.5.25') 
- at g, + sat, 23 ; 25425) 


Eq.4.5.25 suggests a way to reduce the inherent error for 


A@ and af. If g can be estimated in advance, it is better 
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to Solve Eq.4.5.25 than Eq.4.5.10. 
Once we obtaine Aa and Af by Eq.4.5.13 and Eq.4.5.14, 


we have the approximate values of 9, and Ag 


iQ 
N 
i) 


2 eet or = 9, + aa 23. (4.5.26) 
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q nal g, 2.5727) 


oS 
QD 
fl 


Then Eq.4.5.25 can be expressed in the form 


og 


‘5a |, (4.5.28) 


AG At ee aT 


Applying the least-squares procedure to Eq.4.5.28, we can 
get the alternative solution for Af and Aa by using the 
equations which are obtained by replacing g,with 3, in 
meee 5S.13 and &q.4.5.14. The iterative procedure 
calculation of new as and application of LSM - must tend to 


yield more accurate solutions. For example, the iterative 


Af = Af*- aaf 
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Fig24 Ooh snes Schematic diagram showing 
convergence of in iterative computations. 
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calculations for Af are made by 


< 0G.9, >< $433 > - 9) 29. >< 4G, 28 > 


At aa’ da 3a. 
a re ee ek ee et es ae 4.5.29 
’ J Polis | a8 1 ag egy 
£9, 32)<55 , Sr) ee 


ne tendency of .q’- 9 naturally results in af’ aFf*, 
Consequently, Af’ might possibly behave as shown in 
Fig.4.5.1, where g is assumed to be a positive function. 
Another iteration scheme is possible. In Eq.4.5.25', if 
gy,” and Ag’ are respectively substituted for g, and 4g, 
the unknown parameter is only Af. Then the least-squares 
procedure can be applied with respect to Af, and leads to 


the solution 


/ CAG~A F, 93> 
Af- = oe (4.5.30) 
On the other hand, the calculation of A@ is based on the 
equation 
a #23) 
AG z 4 g, + Ad er 5 & tS s ts) 


A@ is calculated by 
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Using Aa’ of Eq.4.5.32, calculate 9/ and 4g’, and repeat 
the computation of Eq.4.5.30 and Eq.4.5.32. The behavior of 


Af’ with the iteration is given by 
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This derives from the substitution of 
RG, PE 9, 7.9 ,95,, 45 =) 9g, 2% (4.5.34) 
mto Eq.4:5.30. 
In order to take account of the higher terms in the 


Taylor's expansion of NG, “Ahq. 2s 5.257, TOLLowing ~ two 


methods are considered. One is based on the equation 





A, (aajy + 2B;)(aaj) + 2C; = 0 (4.5.35) 
where 
ned 3g 
a xe i oe Lj 
ae a5 og (4.5.36) 
fe i era Lj Aleks 
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The above equations are the result of adding the second term 
of the Taylor's expansion to Eq.4.5.25. Since one available 


root cf the above polynomial equation is 


been ( By -— 2A;¢; Ties 


Aaj (4.5.37) 
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The-least squares solution of 4a,’ is given by 
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Rap iape radially, (4.5.38) 
Sills 
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computations, due to mis-estimates in earlier stages of the 
calculations as well as noise disturbance in observed data. 
when Eq.4.5.38 is used, the calculation of 4 in Eq.4.5.36 
should not be made by cee apie. + aa J but by 
a. = 9(a,+48,°)), because Eq.4.5.35 has taken account of a 
higher term. 

Another approach for the treatment of the higher terms 


of Ag is to apply appropriate non-linear inversion methods 


to the equation 


joy) ( . 2 
AG = ati + sat, 23 Psglay Sab e (4.5.39) 


instead of Eq.4.5.28 and/or Eq.4.5.31. However, the usage of 
any non-linear inversion method has no advantage unless it 
can reduce computational effort and yield fast convergence. 
Otherwise, f and a can be calculated by the method described 
in the first part of this section. 

Suppose that two kinds of data are available. There is 
a different way to determine 4@ and Af. Two kinds of data 
will be distinguished by the notation of (1) and (2). We 


have the equation 


(1) ay ag 
AG, Af g da J 
(2) (4.5.40) 
(2) (2) 0S 
" = et AA A 
AG; = 4g; daly 


The solution of this binary simultaneous equation is 
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D; Dj; 


At = 




















Jc ¢gapase teilzee of eetembsee-éim o2 evb .enoisszuqme> | 


yvisedo of sotadsuseisS- spien ec iiev ee enottsluolao | 


‘p oe notratuclad sd2..beeu ef 8€.3.4.p8 nedwy 
a - 


: | See At. & z Pe: yd ebsm ad  3on b{uoda 


a 


st ae SE.2.8.p8 erueced ,(j)1.567,5) 2 4 


mio? t9npim 


ms sipid eds 20 themgesda €f49 TO? A >e07gg6 tenjondé 


bor? sa faisvni teehilmnon @2aL7q°c uggs yiggs o3 2i_ph 


> 


noissups sfig= + 


vad Bojan noleteva: ssenii-non ; 
Sieiy Boe tvolte Ienoizasugmes saves ast 

eb hodism od2 yd Bezeluaiso 4d aso 6 Bae 2 _ seiwieds 
,hel2>ee e280) 70.3264 za zi oid, é 

isliavs $735 #365. 30 abmid ows gad2 s2oqque > 

szeb to ebnid ow? .44 One Sh onimreses 2 yew sneis]aeeam 


sv .{S) Bos |!) 20 acisazom ela yd SedaiupniselS ed lit 


? 


Loy 


where 


(2) (2) 
cr) ag t chi 


Deets —— |, 
4 zs oa ly = oa ly 
bean 
a) 2 
OIp A F209 
Aj Sera AGi xa i (4.5.42) 
{ Cry 2) (tet) 
B, = —( a Se yoy ) 


f 


If appropriate intial values for f and a are specified, Aa 
and Af can be determined by using the observed data at, at 
least, one meaSuring point. If the observed data at several 
measuring points are available, the least-squares solutions 
of 4a and Af are given by 
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(4.5.43) 
Aa = 


The inherent error resulting from Eq.4.5.43 can be reduced 
in the same way as the iterative computations described 


previously. 
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4.6 Observation Equations for Geodetic Measurements 

A horizontal penny-shaped fracture theoretically yields 
a displacement field with radial and vertical components. It 
is a reasonable choice to adopt a geodetic network such that 
measuring points are arranged radially from the well site. 

Let (i) and (j) denote two measuring points located 
radially at distances r, and rj respectively. Also let 
R+ and @, denote the slope distance and the vertical 
angle between these two points at time t from the inception 
of a fracture of radius a and amplitude factor f. The 
displacement vector at the point (i) during the time t will 
be denoted by Wea= f Ure W;)- The change of the 
G@isplacementevector during the-time interval At = t,-.t, 
will be denoted by AU = ( AU; ,4W)). 

First, suppose that we measure the slope distances, R, 


and the vertical angles, 6, at time t=1 and t=2. 





Figv4, St.17e sess Parameters defining 
geodetic measurements and displacement 
vectors at the points (i) and (j). 
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Referring to Fig.4.6.1, we have the relations. 
R2<0s8, ~ R,cos0, = 4Uj} ~ 41h (4.6.1) 
Rasing, - R,sinO, = sWi - awy Ua a2) 

If t=1 is substituted by the time of fracturing start, t=0, 

Eqs.4.6.1 and .2 are given in terms of the absolute 

displacements. 

R,cos6, - R,cos8, = Uj - We ROAR, 
R,sin@, ~R sind, = We - Wi (4.6.4) 
There is no radial displacement at the well site, r=0. 
Therefore,if the point (i) is supposed to be located at r=0, 


Eq.4.6.1 and Eq.4.6.3 are simplified into the form 


Ra ¢os6, — Ri cosh = Aj | (4.6.5) 
R, cos6, — RocosG, = U; (4.6.6) 
When the measurements of R and 9 are given in terms of the 
changes 
Tee Kiel ters y MSR, o>) 4k (4.6.7) 
G5 =e + 26 » I6l>> [46| 
Eq.4.6.1 and Eq.4.6.2 are rewritten as 
aReos0, ~ ABR, sind, = aly - al; (4.6.8) 
aRsin& + 46R,c9s6, = all - awh (4.6.9) 


Let g and h_ be the model functions of the uplift and the 
radial displacement, respectively. We can express the 
displacements and its changes, Ui, Wi,4Ui,4Wi, in the form. 

ay §, h Cae, % > 
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AW, = a5: 9(a4, 7 tia SEY | . 
ass 
SUDStut utd Ons. ofa. Eos 4.6.10: into Eqs.4.6.1 and 2,'Eqs.4.6.3 
and 4, or Eqs.4.6.8 and 9 yields the relations between the 
fracture parameters, (a,f) or (Ada,Af) and the measurements, 
Gh+9 ) or (aR, 40). 


Let 
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From Eqs.4.6.1 and 2, we have 


Rocos®, — Ri cos = Shin, + Aa hi (4.6.12) 

Ra sin, ~ Resin®, =-a$ge - as, 355 (4.6.13) 
From Eqs.4.6.3 and 4, we have 

Rz 2050, - Rocos®, = Syzhis,2 (4.6.14) 

Basin’, — Ko sing, = —4-92) 9 (256. 15) 


Next, suppose that we make measurement of the slope 
distance R alone. Denote the coordinates of two measurement 
points, (i) and (j), prior to the presence of the fracture 
by (r; ,2z,;) and (r,,2)), respectively. Letting 

rij oa Boe) ) Ae 6 212 Fi EST 


S t tS (4.6.16) 
Uz = Up - Uy , Wey WU We 


we have the relations 
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L. apkii = Crt Dy er tee Adis (4.6.18) 


It can be considered practically that 


ee o 2 . 34 


and should be, from the aspect of the theoretical model, 


that 


It follows that, 

bets be tA Us (4.6.19) 
Since the above equation does not involve the vertical 
component of the displacements and makes the problem 
trivial, the following is based on Eq.4.6.18 and the 
equation 


ale be aa ott 


AL=a 
oF 0a 


Then we have an observation equation which relates the 
change of the slope distance 4R and the changes of fracture 


parameters, 4a and 4f, in the form 
e ® 
ReaR= a$CAc his 4 +B 3i5,4) * aa (Ahij + +Bs ij, + pepteca fe. 20%) 


where 


+ na 
Ae = Tig t Ehij (2.6.21) 
BE mst Hy Sy Sig, 
On the other hand, from Eq.4.6.19, we obtain the observation 


equation 


Re4R = rh; 6 aShis 4 d aa Si hiz,t ) peated 
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5. Conclusions 

The application of surface deformation measurement to the 
study of details of induced fracture requires an 
understanding of the accuracy with which surface deformation 
delineates changes of fracture source parameters. These 
include the shape, size, dip, strike of fracture as well as 
the pressure acting on the fracture surface and the 
displacement dislocation across the fracture plane. In order 
to analyse this problem, we need tractable and realistic 
fracture models and appropriate inversion methods. Once 
these tools are available, we can examine surface 
deformation data in terms of the attainable resolution for 
source parameters when noise of defined characteristics is 
superposed on the data. 

Fractures result in several meaSurable quantities such 
as uplift tilt, horizontal displacement and strain. Some of 
these deformation fields are more sensitive to some of the 
source parameters than other fields. The sensitivities, 
which can be calculated from information in chapter 4, will 
indicate in particular cases the optimum observation 
strategy to be used. Particular locations will be sensitive 
to some of the source parameters and particular measurement 
configurations will be sensitive to others. For example a 
Single line of tilt meters distributed in a radial direction 
with a spacing calculated from the results in chapter 4 
would be indicated if equidimensional horizontal fractures 


are expected. The conventional strategy of a circular array 
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does not exploit all possible information. 

A horizontal penny shaped fracture model is not always 
adequate for a study of the form of fractures induced under 
various circumstances. It iS necessary to extend the 
analysis to more general fracture models which have a dip or 
a Shape which is not circular but which are still tractable. 
As long aS a fracture plane can be assumed to be 
characterized by a uniform internal pressure or displacement 
dislocation, the surface deformation field can be expressed 


in the form 
G,sme(n,r; ) 
which consists of a linear parameter, m, depending on the 
internal pressure or the displacement dislocation, and 
several non-linear parameters, tn , depending on geometric 
Source parameters. Even when a fracture does not have a 
plane structure Or uniformity of the pressure or 
displacement dislocation, the surface deformation field can 
be expressed by superposition of the basic model equation, 
in the form 
G, = Lam, 8 My > ry 5 ) 

We can apply an inversion scheme and a method of resolution 
analysis like that described in the previous chapter, to 
these equations. 

The detectability of small magnitude surface 
deformation is the minimum requirement to utilize such 
measurements. Accordingly we have to examine a Survey design 


from the point of view of detectability. Often measuring 







x, ; - 


, , 7 —— 
noitemreink @idiegoc [ls tiolqze Jon agen 
. t 
Se - 
“ = Le DS Bi < 63 1osisor A 


™ 
‘ 
» 
* 
' 
~< 
- 
4 
r 
~ 
; 
_ 
a] 


c As mec 
* a on 
5 7 i 2 . y DSsiz 
- , =~ — 
3 J PBs. 2 f »~norzeaoo 
+ 
- ar 
~< - » aT) 
“ ' 
Tesemgteg 2 o. adetenea od 
—S ‘ 
be al 
-e ) i ait ees © isnisin 
3 
= 5 
7 f —_ r 
4 ‘ T6560 ‘ < at | fy ist 


’ “ Hava .etetemeisg so7vem 


" 
io 4#taaoting 3¢ stusovmse eansig 










a 


yiqqe: neon 


a 





9 bodizem « bae- poteisyni ne 
20 boxe # bre-samien Doierors 


= os ' i —_] « ' Pe 7 fo ; - : 
_ me q =P sy, 5 f ‘Saal D : , a» BP Rt ot 
» > @ Orv 1¢% a9 O80 £3 Ua > 26n7 ~P, 4 : 
~ — | whey view at co 7h oa 2 sl aot y we ees 






164 


points are located at places where the largest ground 
responses are predicted. We may expect that some of the 
source parameters can be estimated by finding the location 
yielding the maximum ( more generally, the extrema or nodes 
) of the deformation field. We can not, however, discuss the 
‘applicability of the deformation measurement from such an 
aspect only. Indeed, a sensitivity analysis shows different 
results for these two cases. In one case, the maximum 
amplitude is compared, and in the other the resolution 
factor is compared. 

Surprisingly, in the literature on hydraulic fracture, 
it is generally assumed that failure will occur by tensile 
fracture. In fact, many in-Situ stress measurements have 
been inferred from hydrofracturing data based on the 
criterion for the occurrence of tensile fracture. There is, 
however, a possibility that fracture takes place as_ shear 
fracture Or aS combination of shear and tensile fractures. 
Basically we have to regard the failure mode as one of the 
unknown source properties. It must be considered that the 
measurement of surface deformation can be used for the 
determination of the failure mode if it is used properly, 
and that the possibility of such an analysis is one of the 
virtues of this geophysical method. 

If we may assume a relatively simple geometry of an 
induced fracture, or if we may infer a rough geometry of a 
fracture in terms of strike and dip, it is mathematically 


possible to estimate the distribution of pressure or 
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displacement dislocation on the fracture plane (or planes), 
by subdividing it (or them) into several grid cells each of 
which can be regarded to _ be ‘Ghat acted ceo by uniform 
pressure or dislocation. 

If the surface deformation is analysed in terms of the 
history of the source parameters, we could obtain 
Significant results which delineate the process of fracture 
growth. We may then expect that these analyses could provide 
important information to substantiate or to revise 
theoretical predictions. Moreover, we expect that the 
surface deformation measurement reveals the in-situ behavior 
of hydraulic fracture which can not be investigated by 
laboratory experiments. For example, the determination of 
in-situ fracture toughness is one of the possible studies. 

Measurements of the terrestrial deformation associated 
with volcanic activity and the occurence of earthquakes can 
also be treated by the techniques developed here. These 
studies differ from the investigation of a man-made 
subsurface deformation due to hydraulic fracturing only in 
scale. The areal extent, magnitude and time scale of 
hydraulic fracture are much smaller than those of earthquake 
and volcanic deformation. | 

Since hydraulic fracturing is conducted at relatively 
shallow depths. it is possible to obtain information about 
fracture parameters by direct measurements using observation 
wells. This implies that we can use more control data for 


the analysis of measured surface deformations. Results from 
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the surface deformation analysis supported by such an 
advantage must play an important role not only in the field 
of energy resource development but also in geophysics 


itself. 
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